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Abstract 

Several  observations  of  atmospheric  turbulence  statistics  have  been  reported 
which  do  not  obey  Kolmogorov’s  power  spectral  density  model.  These  observations 
have  prompted  the  study  of  optical  propagation  through  turbulence  described  by 
non-classical  power  spectra.  This  thesis  presents  an  analysis  of  optical  propaga¬ 
tion  through  turbulence  which  causes  index  of  refraction  fluctuations  to  have  spatial 
power  spectra  that  obey  arbitrary  power  laws.  The  spherical  and  plane  wave  struc¬ 
ture  functions  are  derived  using  Mellin  transform  techniques  and  are  applied  to  the 
fleld  mutual  coherence  function  (MCE)  using  the  extended  Huygens-Fresnel  princi¬ 
ple.  The  MCF  is  used  to  compute  the  Strehl  ratio  of  a  focused,  constant  amplitude 
beam  propagating  in  non-Kolmogorov  turbulence  as  the  power  law  for  the  spectrum 
of  the  index  of  refraction  fluctuations  is  varied  from  -3  to  -4.  The  relative  contri¬ 
butions  of  the  log  amplitude  and  phase  structure  functions  to  the  wave  structure 
function  are  computed.  If  inner  and  outer  scale  effects  are  neglected,  no  turbulence 
exists  when  the  power  law  equals  -3.  At  power  laws  close  to  -3,  the  magnitude  of  the 
log  amplitude  and  phase  perturbations  are  determined  by  the  system  Fresnel  ratio. 
At  power  laws  approaching  -4,  phase  effects  dominate  in  the  form  of  random  tilts. 
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LASER  BEAM  PROPAGATION 
IN 


NON-KOLMOGOROV 
ATMOSPHERIC  TURBULENCE 

/.  Introduction 

1.1  Background 

It  is  well  known  that  atmospheric  turbulence  limits  the  performance  of  imaging 
and  laser  systems.  Temperature  variations  associated  with  turbulent  eddies  cause 
random  variations  in  the  index  of  refraction.  This  random  variation  in  the  index 
of  refraction  distorts  an  optical  wave  as  it  propagates  through  the  turbulent  atmo¬ 
sphere.  The  distortions  in  the  optical  wave  caused  by  the  random  index  of  refraction 
lead  to  significant  blurring  in  imaging  systems  and  increased  scintillation  (twinkling) 
and  wander  in  a  laser  beam. 

To  date,  the  performance  of  imaging  and  laser  systems  has  been  estimated  as¬ 
suming  the  Kolmogorov  model  for  atmospheric  turbulence.  The  Kolmogorov  model 
is  a  statistical  model  for  the  variations  of  the  index  of  refraction  in  the  atmosphere. 
The  model  is  defined  by  a  3-dimensional  power  spectral  density.  While  the  Kol¬ 
mogorov  model  has  shown  good  agreement  with  experiment  in  the  past,  recent  ev¬ 
idence  has  shown  significant  deviations  from  the  Kolmogorov  model  in  certain  por¬ 
tions  of  the  atmosphere  [1,  4,  6,  7].  The  power  spectrum  of  turbulence  in  portions  of 
the  troposphere  and  stratosphere  may  exhibit  non- Kolmogorov  statistics  [6,  28].  The 
propagation  of  light  through  turbulence  which  does  not  obey  Kolmogorov  statistics 
is  not  well  understood. 
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1.1.1  Air  Force  interests  in  non- Kolmogorov  turbulence.  Understanding 
the  propagation  of  light  through  non-Kolmogorov  turbulence  is  important  to  two 
distinct  USAF  missions: 

1.  Theater  ballistic  missile  defense. 

Phillips  Laboratory  has  proposed  developing  a  prototype  airborne  laser  (ABL) 
system  capable  of  destroying  a  SCUD-type  missile  at  very  long  ranges.  The 
system  is  intended  to  operate  in  the  troposphere  and  stratosphere,  and  requires 
laser  energy  to  be  transmitted  to  the  target  over  long,  horizontal  paths.  This 
is  precisely  where  non- Kolmogorov  turbulence  has  been  observed  [28].  Under¬ 
standing  the  effects  of  non- Kolmogorov  turbulence  on  laser  beam  propagation 
may  play  a  key  role  in  determining  the  operational  limits  of  the  ABL  concept. 

2.  Ground-based  imaging  of  space  satellites. 

Air  Force  Space  Command  has  a  mission  to  gather  intelligence  on  foreign  satel¬ 
lites  in  orbit.  The  mission  is  accomplished  via  the  use  of  ground  based  tele¬ 
scopes.  Unfortunately,  turbulence  limits  the  resolution  that  a  telescope  can 
achieve.  Some  turbulence  effects  can  be  overcome  through  the  use  of  adaptive 
optics.  The  performance  of  adaptive  optical  systems  can  be  further  improved 
by  using  prior  knowledge  of  turbulence  statistics  (e.g.,  Kolmogorov  statistics). 
If  the  turbulence  is  modeled  as  Kolmogorov  and  is  actually  non- Kolmogorov 
through  a  portion  of  the  atmosphere,  the  performance  of  these  adaptive  optical 
systems  will  be  degraded. 

The  research  presented  in  this  thesis  is  primarily  devoted  to  addressing  the  problems 
stated  above,  but  has  potential  application  in  the  areas  of  laser  radar  systems  and 
other  remote  sensing  applications,  airborne  and  earth-to-space  laser  communications 
systems,  and  satellite  power  beaming  systems. 

1.1.2  The  nature  of  refractive  turbulence  and  its  power  spectrum.  Refrac¬ 
tive  turbulence  is  generated  by  the  differential  heating  of  the  Earth’s  atmosphere. 
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winds,  and  large  scale  weather  phenomena.  As  large  pockets  (or  eddies)  of  warm  air 
mix  with  colder  air,  energy  is  transferred  from  the  larger  eddies  to  eddies  of  smaller 
sizes.  These  smaller  eddies  transfer  their  energy  to  still  smaller  eddies.  This  en¬ 
ergy  cascade  continues  until,  at  some  small  scale,  viscous  effects  dominate,  the  flow 
becomes  laminar  and  energy  is  dissipated  as  heat. 

In  the  seminal  work  on  the  statistics  of  atmospheric  turbulence,  Kolmogorov 
[15]  derived  the  structure  function  for  the  velocity  of  a  turbulent  flow, 

D,(r)  =  {[Kri  +  r)  -  t.(r01")  =  Cy!^  L<r<  L„  (1.1) 

where  the  angle  brackets  denote  the  statistical  expectation  operator,  'y(ri)  is  the 
velocity  vector  at  point  ri,  is  the  velocity  structure  constant,  the  scalar  r  is  the 
magnitude  of  the  vector  r,  and  and  Lg  are  the  inner  and  outer  scales,  respectively. 
The  domain  of  r  between  Ig  and  Lg  is  known  as  the  inertial  subrange  (See  Fig. 
1.1).  Tatarski  [25]  used  the  concept  of  a  conservative  passive  additive  to  relate  the 
velocity  structure  function  to  the  structure  function  for  potential  temperature  and 
then,  using  the  Gladstone-Dale  relation  [5,  21],  derived  the  structure  function  for 
the  variations  in  the  index  of  refraction, 

Dn(r}  =  {|n(ri  +  r)  -  n(r,)|")  =  Cy/^  !.  <  r  <  (1.2) 

where  n(ri)  is  the  index  of  refraction  at  the  point  ri  and  is  the  index  structure 
constant.  It  is  shown  in  Chapter  2  that  Eq.  (1.2)  implies  the  random  variations  in 
the  index  of  refraction  obey  a  three  dimensional  power  spectrum  given  by 

$„(k,z)  =  0.033C2(z)k-ii/3  27r/Lo  <  k  <  27r/lo,  (1.3) 

where  $„(«,  z)  is  the  power  spectral  density  as  a  function  of  position  along  the  optical 
path  z,  K  is  the  spatial  wavenumber  and  Cl{z)  is  the  index  structure  constant  which 
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is  allowed  to  vary  smoothly  along  the  optical  path.  C^{z)  describes  the  strength  of 
the  turbulence  and  has  units  of  Equation  (1.3)  is  known  as  the  Kolmogorov 

power  spectrum.  The  Kolmogorov  spectrum  agrees  well  with  experiment  for  vertical 
propagation  paths  [21]. 

The  key  assumptions  in  Kolmogorov’s  and  Tatarski’s  derivations  are  those 
of  homogeneity  and  isotropy.  For  our  purposes,  “homogeneous”  turbulence  means 
the  turbulence’s  power  spectrum  does  not  depend  on  location  in  the  atmosphere. 
“Isotropic”  turbulence  means  the  turbulence’s  power  spectrum  does  not  depend  on 
the  direction  in  which  the  optical  wave  is  travelling.  The  Kolmogorov  power  spec¬ 
trum  is  valid  in  the  inertial  subrange  (see  Fig.  1.1),  but  clearly  breaks  down  for 
physical  reasons  at  the  inner  and  outer  scales.  At  the  inner  scale  (smallest  eddy 
sizes)  the  turbulent  flow  becomes  laminar  and  the  kinetic  energy  of  the  flow  is  dis¬ 
sipated  as  heat.  At  the  outer  scale  (largest  eddy  sizes)  the  spectrum  should  be 
limited  by  the  flnite  size  of  the  atmosphere,  thus  the  assumptions  of  homogeneity 
and  isotropy  are  violated.  Fortunately,  it  has  been  shown  that  many  calculations 
of  interest  in  propagation  theory  are  insensitive  to  the  exact  nature  of  the  power 
spectrum  at  the  inner  and  outer  scales. 


1.2  The  Existence  of  Non- Kolmogorov  Turbulence 

Von  Karman  [27],  upon  examination  of  Kolmogorov’s  assumptions,  proposed 
an  alternative  power  spectrum  which  addresses  the  problems  at  the  inner  and  outer 
scales,  but  retains  the  11/3  power  law  in  the  inertial  subrange.  Von  Karman  proposed 
the  following  power  spectrum. 


#„(«,z)  =  0.033C2(.e:)(«2  + 


Ko  <  <  K'n 


(1.4) 


where  Kq  =  27r/Lo  and  =  2t: flo-  In  many  optical  calculations  involving  turbu¬ 
lence,  the  von  Karman  spectrum  is  used  for  mathematical  convenience  because  it 
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does  not  diverge  to  infinity  at  zero  wavenumber  (infinite  eddy  size).  Equation  (1.4) 
also  falls  off  rapidly  at  the  inner  scale.  The  attribute  of  finite  energy  coupled  with 
a  rapid  decrease  in  energy  for  increasing  wavenumber  allows  many  of  the  integrals 
involved  in  turbulence  theory  to  be  evaluated  numerically.  Figure  1.1  shows  the 
differences  between  the  Kolmogorov  and  von  Karman  power  spectra. 


Kolmogorov  Spectrum 
Von  Karman  Spectrum 


Figure  1.1  Kolmogorov  and  von  Karman  power  spectra  for  =  10“^^m“^/^.  The 
plot  is  in  terms  of  spatial  wavenumber  or  27r/eddy  size,  thus  the  inner 
scale  (smallest  eddies)  are  to  the  right  and  the  outer  scale  (largest  ed¬ 
dies)  are  to  the  left.  The  inertial  subrange  is  the  constant  slope  portion 
of  the  von  Karman  spectrum. 


Other  authors  have  proposed  extensions  of  Kolmogorov’s  work.  Greenwood 
[14]  and  Hill  [13]  have  found  a  slight  peak  before  the  exponential  decay  of  the  inner 
scale.  Gurvich  [11]  has  proposed  a  local  axisymmetric  model  in  place  of  the  isotropic 
model  to  account  for  the  existence  of  buoyancy  forces  due  to  temperature  differentials 
between  the  turbulent  eddies. 
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In  addition  to  these  theoretical  works,  several  experimental  observations  have 
exhibited  peculiar  power  spectra.  Buser  [4]  has  made  measurements  of  turbulence 
power  spectra  in  the  Earth’s  boundary  layer  by  using  a  3-beam  Mach-Zender  opti¬ 
cal  interferometer.  By  interfering  several  laser  beams  propagating  at  different  paths 
through  the  turbulence  and  analyzing  the  interferograms,  Buser  found  that  the  Kol¬ 
mogorov  theory  did  not  support  his  measurements.  The  statistics  still  obeyed  a 
power  law  but  the  observed  exponent  was  significantly  lower  than  11/3.  Buser  iden¬ 
tified  2  possible  explanations  for  this: 

1.  At  ground  level  where  these  measurements  were  conducted,  the  conditions  for 
vertical  homogeneity  and  isotropy  are  most  probably  violated,  especially  at 
larger  scale  sizes. 

2.  The  buoyancy  forces  mentioned  earlier  may  be  large  enough  to  effect  the  results 
of  the  experiment. 

Dayton,  et.  al.  [6],  have  measured  non- Kolmogorov  phase  structure  functions 
with  a  Shack-Hartmann  wavefront  sensor.  These  measurements  imply  the  existence 
of  non- Kolmogorov  turbulence  power  spectra. 

Dalaudier  [7]  has  found  direct  evidence  of  temperature  sheets  in  the  upper 
atmosphere.  These  sheets  are  areas  of  strong  laminar  flow  having  significant  tem¬ 
perature  differentials  with  the  surrounding  air,  thus  the  conditions  of  isotropy  and 
homogeneity  are  not  necessarily  valid  in  the  upper  atmosphere. 

Additionally,  the  existence  of  non-Kolmogorov  turbulence  is  supported  by 
Wissler’s  [28]  preliminary  analysis  of  the  ARGUS  anemometer  data.  The  ARGUS 
anemometer  is  a  high  frequency  (12kHz)  temperature  probe  mounted  near  the  stag¬ 
nation  point  on  the  nose  of  a  modified  RC-135  aircraft.  ARGUS  has  been  tasked 
by  Phillips  Laboratory  to  measure  refractive  turbulence  spectra  at  high  altitudes  for 
the  Airborne  Laser  Program  (ABL).  The  index  of  refraction  at  high  altitude  can  be 
estimated  by  measuring  temperature  and  pressure  as  a  function  of  the  flight  path 
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and  relating  the  measured  temperature  and  pressure  to  refractive  index  through  the 
Gladstone-Dale  relationship.  Many  observations  of  non- Kolmogorov  turbulence  are 
believed  to  exist  in  the  ARGUS  data.  However,  due  to  the  difficulties  associated  with 
the  frequency  response  calibration  of  the  anemometer  probes,  the  ARGUS  data  is 
inconclusive  at  this  time. 

Because  of  the  increasing  number  of  experimental  observations  of  non- Kolmogorov 
turbulence,  it  has  become  necessary  to  develop  a  theory  which  provides  a  basis  for 
estimating  the  performance  of  optical  systems  in  non- Kolmogorov  turbulence.  For 
purposes  of  this  thesis,  non- Kolmogorov  turbulence  is  defined  as  turbulence  whose 
three  dimensional  power  spectrum  obeys  an  arbitrary  power  law, 

a,  z)  =  a(a)/?(2)«~“,  (1.5) 

where  a,  z)  is  the  power  spectral  density  as  a  function  of  position  along  the 
optical  path  2,  k  is  the  spatial  wavenumber,  a  is  the  power  law,  jd(z)  is  the  index 
structure  constant  along  the  path  of  propagation  and  has  units  of  m^~“,  and  a(a) 
is  a  function  which  maintains  consistency  between  the  index  structure  function  and 
its  power  spectrum  (See  Eq.  2.49). 

1.3  Methods  for  Analyzing  the  Ejfects  of  Turbulence  on  Laser  Beam  Propagation 

In  order  to  determine  a  suitable  method  to  analyze  the  effects  of  atmospheric 
turbulence,  one  must  establish  criteria  for  suitability.  For  purposes  of  this  research, 

I  use  the  following  criteria; 

1.  The  chosen  method  must  accommodate  solutions  for  an  arbitrary  wave,  not 
just  a  plane  or  spherical  wave. 

2.  The  chosen  method  must  accommodate  solutions  for  a  turbulence  power  spec¬ 
trum  which  obeys  an  arbitrary  power  law. 
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3.  The  chosen  method  must  provide  good  physical  insight  to  the  problem.  That  is, 
the  solution  method  must  separate  the  geometry  of  the  problem,  the  statistics 
of  the  atmosphere,  and  the  characteristics  of  the  optical  wave  to  be  propagated. 

Lee  and  Harp  [16],  solved  the  plane  wave  and  spherical  wave  cases  by  breaking 
up  the  turbulent  atmosphere  into  a  sum  of  thin  slabs  perpendicular  to  the  propaga¬ 
tion  path.  Each  thin  slab  has  its  two  dimensional  refractivity  field  decomposed  into 
its  Fourier  components.  The  effect  of  each  Fourier  component  of  the  turbulence  acts 
as  a  phase  diffraction  grating  on  the  plane  wave.  Once  these  effects  have  been  com¬ 
puted,  the  authors  sum  the  contributions  from  each  Fourier  component  in  each  slab 
statistically  over  the  path.  This  approach  yields  a  solution  of  three  multiplicative 
terms: 

1.  A  term  containing  the  power  spectrum  of  the  turbulence. 

2.  A  term  relating  the  fluctuations  of  a  wave  at  a  point  in  the  receiver  plane  with 
respect  to  another  point  in  the  receiver. 

3.  A  term  which  relates  the  relative  efficiency  of  a  turbulent  eddy  of  a  certain 
size,  at  a  certain  point  in  the  path,  in  producing  perturbations  in  the  receiver 
plane.  This  term  is  known  as  a  “filter”  function. 

This  method  yields  more  physical  insight  to  the  problem,  but  is  still  limited  to  the 
plane  and  spherical  wave  cases. 

Lutomirski,  Yura,  and  Buser  [17,  18,  30],  derived  a  method  for  propagating 
arbitrary  fields  through  turbulence,  using  the  definition  of  the  mutual  coherence 
function  (MCF)  and  the  extended  Huygens-Fresnel  principle.  The  mutual  coherence 
function  of  a  random  wave  is  the  ensemble  average  of  the  complex  amplitude  of  the 
electric  field  at  one  point  multiplied  by  the  complex  amplitude  of  the  field  at  another 
point.  Mathematically  the  MCF  is  given  by, 

MCF{xi,yi,X2,y2)  =  {U{xi,yi)U*{x2,y2)),  (1-6) 
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Figure  1.2  A  wave  propagating  from  the  ^,77  plane  to  the  x,y  plane,  using  the 
extended  Huygens-Fresnel  principle.  The  wavefront  is  the  summation 
of  the  envelopes  of  spherical  wavelets  whose  centers  where  located  on 
an  earlier  wavefront. 


where  U{xi,yi)  and  U{x2,y2)  are  the  complex  amplitudes  of  the  field  at  the  points 
{xi,yi)  and  (x2,y2)  and  the  angle  brackets  denote  the  ensemble  average.  The 
Huygens-Fresnel  principle  states  (see  Fig.  1.2): 

Each  point  on  an  optical  wave  can  be  regarded  as  the  center  of  a  sec¬ 
ondary  disturbance  which  gives  rise  to  spherical  wavelets  and  the  position 
of  the  wavefront  at  any  later  time  is  the  envelope  of  all  such  wavelets. 

Thus,  if  one  knows  how  turbulence  effects  a  spherical  wave,  one  can  determine  the 
effects  of  turbulence  on  an  arbitrary  wave.  The  MCF  approach  provides  a  solution 
to  the  general  problem  and  provides  a  clear  separation  between  the  geometry  of  the 
problem,  the  statistics  of  the  atmosphere,  and  the  wave  to  be  propagated. 

Beland  [1],  in  an  unpublished  work,  has  used  the  methods  and  results  of 
Tatarski  [25]  and  Sasiela  [23]  to  analyze  the  incoherent  imaging  case  where  the 
turbulence  spectrum  is  non-Kolmogorov.  In  this  work,  Beland  derives  the  general 
expression  for  the  wave  structure  function  for  atmospheric  turbulence  which  obeys 
an  arbitrary  power  law. 

By  combining  Beland’s  structure  function  calculation  [1]  with  the  methodology 
of  Lutomirski  and  Buser  [17],  a  method  is  developed  for  propagating  an  arbitrary 
wave  in  atmospheric  turbulence  which  obeys  an  arbitrary  power  law. 
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1.4  Organization 

In  Chapter  2,  the  mathematical  background  required  for  this  thesis  is  intro¬ 
duced.  Mellin  transform  techniques  are  briefly  reviewed.  Next,  the  relationship 
between  the  index  power  spectrum  of  turbulence  and  the  index  structure  function 
is  examined.  Next,  the  spherical  wave  structure  function  for  non- Kolmogorov  tur¬ 
bulence  is  derived  and  a  generalized  expression  for  Fried’s  coherence  diameter  To  is 
found.  Finally,  an  expression  for  the  variance  of  the  log  amplitude  fluctuations  of  an 
optical  wave  is  defined  and  the  limitations  of  weak  fluctuation  theory  are  examined 
for  various  power  laws. 

Chapter  3  describes  the  extended  Huygens- Fresnel  principle  for  propagation 
of  light  in  a  vacuum  and  in  a  random  medium.  The  mutual  coherence  function  of 
an  arbitrary  wave  in  a  random  medium  is  derived.  The  MCF  is  then  used  to  calcu¬ 
late  the  on-axis  intensities  of  an  arbitrary  wave  in  vacuum  and  in  non- Kolmogorov 
turbulence. 

In  Chapter  4,  an  expression  for  the  Strehl  ratio  of  a  focused,  constant  amplitude 
wave  is  derived.  The  asymptotic  behavior  of  Fried’s  resolution  metric  7^/7^max  is 
examined  for  various  power  laws.  Additionally,  several  comparisons  of  Strehl  ratio 
as  a  function  of  the  power  law  are  made  under  various  constraints. 

In  Chapter  5,  the  relative  contributions  of  the  log  amplitude  and  phase  struc¬ 
ture  functions  to  the  wave  structure  function  are  addressed. 

Chapter  6  discusses  the  major  findings  of  this  research  and  recommendations 
for  future  research. 
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II.  Mellin  Transforms  and  the  Wave  Structure  Function 

In  order  to  understand  the  mathematical  development  of  the  non- Kolmogorov 
propagation  theory,  one  must  be  familiar  with  Mellin  transform  techniques.  The 
Mellin  transform  is  a  mathematical  tool  which  allows  many  complicated  real  and 
complex  integrals  to  be  evaluated  in  a  straightforward  and  efficient  manner.  The 
first  section  of  this  chapter  is  dedicated  to  a  short  review  of  Mellin  transforms.  After 
the  section  on  Mellin  transforms,  the  relationship  between  the  index  power  spectrum 

O',  2)  and  the  index  structure  function  T>„(r)  is  derived.  This  relationship 
allows  one  to  understand  why  the  power  law  a  of  the  index  power  spectrum  can  only 
take  on  certain  values.  In  Section  2.3,  the  optical  wave  structure  function  is 

calculated  from  the  spherical  wave  correlation  function  which  can  be  derived  from 
the  method  of  small  perturbations.  The  wave  structure  function  will  be 

used  extensively  in  Chapters  3  and  4  to  calculate  the  mutual  coherence  function 
(MCF)  and  Strehl  ratio  of  a  laser  beam.  Since  many  of  the  results  contained  in  this 
research  are  normalized  to  atmospheric  and  path  conditions,  a  generalized  version 
of  Fried’s  coherence  diameter  is  derived.  Finally,  the  limitations  of  the  non- 
Kolmogorov  theory,  imposed  by  the  method  of  small  perturbations  [2,  5,  25]  are 
examined.  The  limitations  are  examined  by  constraining  the  variance  of  the  log 
amplitude  fluctuations  of  the  optical  field  to  the  regime  in  which  the  method  of 
small  perturbations  is  known  to  agree  with  experiment. 

2.1  Introduction  to  Mellin  Transforms 

In  many  problems  in  propagation  theory  we  are  concerned  with  integrals  of 
the  form 


h{a)  =  d«K$„(K,a,z)/i(oi/<;*^)/2(a2K^").../n(a™K''"),  (2.1) 


2-1 


where  a,  z)  is  the  3-dimensional  power  spectrum  of  the  turbulence  and  ) 

are  functions  of  the  integration  variable.  Because  the  Mellin  transform  defined  in 
Eq.  (2.2)  has  a  form  similar  to  the  integral  above  (Eq.  (2.1)),  the  Mellin  transform 
is  a  very  useful  tool  for  evaluating  integrals  associated  with  propagation  of  light 
through  turbulence.  In  this  section,  the  Mellin  transform  is  defined  and  some  of 
its  properties  are  proven.  Euler’s  gamma  function  is  discussed  and  some  elemen¬ 
tary  Mellin  transforms  are  evaluated  as  examples.  For  more  detailed  discussions  on 
Mellin  transforms,  one  should  refer  to  the  works  of  Sasiela  [22,  23,  24]  and  Marichev 
[19]. 

The  Mellin  transform  and  its  inverse  are  defined  as 

roo 

H  (s)  =  M[h{x)]  =  /  dxx^'~^h{x)  (2.2) 

J  0 

and 

h{x)  =  M-^[H(s)]  =  (2.3) 

where  the  path  of  integration  for  the  inverse  transform  is  a  straight  line  from  -joo 
to  joo.  The  contour  c  crosses  the  real  axis  at  a  value  of  s  for  which  the  forward 
transform  converges.  This  strip  of  the  complex  plane  for  which  the  contour  c  crosses 
the  real  axis  is  the  domain  of  the  transform  V. 

2.1.1  Mellin  transform  properties.  The  Mellin  transform  possesses  several 
properties  which  can  be  used  to  extend  a  table  of  Mellin  transforms.  In  this  section, 
the  following  Mellin  transform  properties  are  proven. 


M[h{ax)]  - 

a  ^H{s),  a  >  0 

(2.4) 

M[x'^h{x)]  - 

i?(s  -|-  g) 

(2.5) 

M[h{x'^)] 

(2.6) 
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Theorem  II. 1  Let  h{x)  he  a  Mellin  transformable  function  on  domain  V,  then 


M.[h{ax)]  — >  a  ^H{s)  a  >  0. 

(2.7) 

Proof: 

/•OO 

M[h{ax)\—  /  dxx  h{ax) 

JO 

(2.8) 

Let 

r  —  ax  dr  =  adx 

X  =  rja  da;  =  dr/a 

then, 

(2.9) 

Thus, 

roo 

M[h{ax)]  =  a~^  /  drr"*  ^h{r) 

J  0 

(2.10) 

and 

fA\h{ax)\  — ^  a~^H{s)  a  >  0. 

(2.11) 

Theorem  11.2  Let  h{x)  be  a  Mellin  transformable  function  on  domain  V,  then 


M.[x°'h{x)]  H{s  +  a). 


(2.12) 


Proof: 


fOO 

M[x‘'h{x)]  —  j  da:  a;*  ^x°'h{x) 


(2.13) 


Let  r  =  s  +  a,  then 


fOO 

M{x°'h{x)]  =  /  dx  x'^~^h{x)  =  H{r), 

J  0 


(2.14) 
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thus, 


M[x°'h{x)]  —^H{s  +  a). 


(2.15) 


Theorem  II.3  Let  h{x)  he  a  Mellin  transformable  function  on  domain  V,  then 


M[h{x^)]  - 

P#0. 

\p\ 

(2.16) 

Proof: 

M[h{x^)]  ■■ 

/‘OO 

=  /  da;  a;*  ^h{x^) 

(2.17) 

Let 

r  =  x^ 

dr  =  pa;^“^da; 

X  — 

dr  — 

then, 

M[h{x^)]  =  j 

i  '  M--). 

(2.18) 

Thus, 

M[h{x^)]  - 

1  a 

=  -  dr  rp~^h{r) 

p  Jo 

(2.19) 

and 

M[h{x‘^)]  - 

m 

(2.20) 

2.1.2  Gamma  functions  and  Mellin  transform  examples.  Mellin  transforms 
can  typically  be  expressed  as  ratios  of  Euler’s  gamma  functions.  The  gamma  function 
is  defined  as 

roo  °°  t  — 1  r°° 

r{s)=[  da;exp(-a:)a:*“^  =  £  p - T“  +  /  exp(-a:)a;^  .  (2.21) 

Jo  n=0  s  -r  n  Ji 
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The  last  integral  is  analytic  on  the  entire  complex  plane.  The  gamma  function 
possesses  simple  poles  at  the  negative  integers  with  residues  (-l)”/n!.  To  express 
ratios  of  gamma  functions  the  notation  of  Sasiela  [23]  is  used, 

^  QI^  0-2  •  *  • 

/?!  (^2  ...  Pn 

Most  of  the  integrals  in  this  thesis  can  be  evaluated  by  a  simple  change  of  variables 
and  a  table  lookup,  just  as  in  Fourier  transform  theory.  Tables  can  be  found  in  Sasiela 
[22,  23,  24]  and  Marichev  [19].  The  example  below  is  from  a  problem  discussed  in 
the  following  section. 


r[c>!i]r[Q;2] . . .  r[Q!m] 

mim] . . .  r[/?™]  • 


Example  II.  1  In  the  Section  2.2,  we  wish  to  find  the  power  spectrum  associated 
with  a  structure  function  that  obeys  an  arbitrary  power  law.  Tatar  ski  derived  the 
following  relation, 


^n(«) 


1 


sin(Kr)  d 
Kr  dr 


dr. 


(2.23) 


where  *5  the  isotropic  3-dimensional  turbulence  power  spectrum  and  Dn{r) 

is  the  index  structure  function.  Substituting  the  index  structure  function  D„(r)  = 
I3{z)r'^  (See  Eq.  (2.52))  into  Eq.  (2.23)  and  evaluating  the  derivatives  yields 


^n{K,z) 


/?(z)(7"  +  7) 

Jo 


Kr 


(2.24) 


Making  a  change  of  variables 


X  —  Kr  da;  =  «dr 
r  =  x/k  dr  =  dx/K, 


2-5 


Eq.  (2.24)  becomes 


/?(2:)(7^  +  7)  r 

47r2^3+7  Jq 


da:  a;'*'  ^ 


sin(a:). 


(2.25) 


By  defining  a  =  'j  +  3,  Eq.  (2.27)  becomes, 

P 

z)  =  2“-®/?(2)(a^  -  5a  +  6)7r~^/^  2  <  a  <  4.  (2.28) 

Notice  the  domain  in  Eq.  (2.28)  differs  from  Eq.  (2.54)  ^<1-  (2.57).  Equation 

(2.28)  must  be  analytically  continued  to  the  domain  3  <  a  <  5,  to  ensure  the  power 
spectral  density  is  positive.  The  analytic  continuation  is  done  automatically  in  the 
next  example. 

In  many  turbulence  problems  it  is  necessary  to  evaluate  the  Mellin  transform  of  a 
function  minus  the  first  term  of  its  power  series.  The  resultant  transform  is  the  same 
as  the  original  function  but  the  domain  of  the  transform  has  moved  past  one  pole. 
This  is  demonstrated  in  the  following  example. 
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Example  II. 2  To  derive  Beland’s  [1]  expression  for  a{a)  in  Eq.  (1.5)  we  can  use 
Eq.  (2.48)  which  is  derived  later  in  the  chapter, 


Dn{r,z)  =  Stt  y  d/c 


1  - 


sin(«;r) 


Kr 


(2.29) 


Substituting  the  definitions  for  the  index  structure  function  Z)„(r,  z)  =  P{z)r'^  and 
turbulence  power  spectrum  2)  ==  a{a)p{z)K,~°‘  (See  Eq.  (2.52)  and  Eq.  (2.56)) 
into  Eq.  (2.29)  yields 


—  87ra(Q:)  J 


sin(/cr) 

/  dKK^ 

1  ^  ' 

Jo 

Kr 

Letting  s  =  3  —  ck  and  making  the  change  of  variables, 

X  =  Kr  ds  =  rd« 
n  =  x/r  dK  =  dx/r, 


roo 

/  dx  x^  ^ 

- 1 

1 

_ 1 

Jo 

L  ®  J 

Eq.  (2.30)  becomes 

=  8TTa(a)r°‘~' 

Multiplying  Eq.  (2.32)  by —1  yields 

fOO  r 

r'^  =  — 87ra(a)r““^  /  dx  x”"^  x~^  sin(a;)  —  1 
Jo  '• 

Applying  the  Mellin  transforms 


M[x°'h{x)] 

Ad[sin(a;)]  ^  2*“^V^r 


H{s  A  a) 

1/2  + s/2 

I -s/2 


|^(s)|  <  1 


(2.30) 


(2.31) 


(2.32) 


(2.33) 
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M\l]  ->  lim,^o  1S(«)I  <  « 


(2.34) 


yields 


— 87ra(a)r" 

(3-s)/2 


1  1 

- 1 - 

s  +  e  s  —  e 


|3fJ(s)|  <  e.  (2.35) 


In  the  limit  as  e  0  the  poles  at  0  and  -e  additively  cancel  and  the  domain  becomes 
-2  <  3?(s)  <  0.  Since  s  =  S  -  a,  Eg.  (2.35)  becomes 


,  ,  ^  (3-a)/2 

r'*  —  —8'jra{a)r°'  ^2^  3  <  a  <  5. 

a  12 


(2.36) 


Equating  the  exponents  of  r  and  solving  for  a{a)  yields 


a(a)  =  -(2)“-"7r-3/2r  '  3  <  a  <  5. 

(3  -  a)/2 


(2.37) 


The  moving  of  a  domain  past  a  pole  is  an  example  of  an  analytic  continuation  of  a 
function. 


2.2  The  Relationship  Between  the  Index  Power  Spectrum  and  the  Index  Structure 
Function 

The  random  variations  in  the  index  of  refraction  can  be  described  by  the  three 
dimensional  autocorrelation  function 


r„(rl,rl)  =  £  {n(rl)n(rl)}  , 


(2.38) 
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where  n(rl)  and  n(rl)  are  the  indices  of  refraction  at  points  r1  and  r"^  and  £  is  the 
expectation  operator.  By  assuming  that  the  index  variations  are  homogeneous  (sta¬ 
tionary)  throughout  the  atmosphere,  the  autocorrelation  function  can  be  expressed 
as  a  function  of  difference  coordinates  only, 

r„(r)  =  £  {n(r;)n(rl  -  r)}  ,  (2.39) 

where  r  =  -  rl.  The  Weiner- Khinchine  theorem  states  the  autocorrelation  and 

power  spectral  density  of  a  wide  sense  stationary  random  process  are  Fourier  trans¬ 
form  pairs,  therefore 

Assuming  the  turbulence  is  isotropic  (i.e.,  r„(f)  depends  only  on  r  =  |r|),  Eq.  (2.40) 
can  be  written  as  a  single  integral  [3], 

1 

$„(«)  =  rTn{r)sin{Kr)dr  (2.41) 

and  the  inverse  relationship  is 

Air 

r,i(r)  =  —  /  k$„(k)  sin(K7-)dK.  (2.42) 

r  Jq 

For  isotropic  turbulence,  Eq.  (2.41)  and  Eq.  (2.42)  are  the  fundamental  rela¬ 
tionship  between  the  autocorrelation  of  index  variations  and  the  index  power  spec¬ 
trum.  It  is  through  the  use  of  Eq.  (2.41)  that  the  index  structure  function,  defined 
below,  is  related  to  the  index  power  spectrum. 

An  extremely  useful  method  of  describing  the  random  variations  in  the  index 
of  refraction  is  through  the  use  of  structure  functions.  The  structure  function  of  the 
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index  of  refraction  variations  is  defined  as 


^n(rl,rl)  =  5  (Krl)  -  ^(ra)]^}  .  (2.43) 

Expanding  Eq.  (2.43)  and  applying  the  linearity  of  the  expectation  operator  yields 
Dn{ri,f*2)  =  S  {n(rl)2}  -  28  {n(r;)n(r2)}  +  8  {n(r^)2}  .  (2.44) 

Since  the  turbnlence  is  assumed  wide  sense  stationary, 

e  {n(rl)"}  =  S  {n(rlf }  =  r.{0).  (2.45) 

Thus,  the  index  structure  function  for  isotropic  turbulence  can  be  written 

D„(r)  =  2[r,(0)  -  r„(r)].  (2.46) 


Substituting  Eq.  (2.42)  into  Eq.  (2.46)  yields 


r  ( 47r  /“^  1  47r 

Dn{r)  =  2  lim| — J  «$„(«;)  sin(«;r)dK| «$„(«)  sin(Kr)dK  ,  (2.47) 


where  the  limit  on  the  first  term  is  required  because  of  the  pole  at  r  =  0.  Evaluating 
the  limit  yields 


Dn{r)  =  Stt  /  $„(«;)«; 

V  0 


sin(«r) 


Kr 


dK. 


(2.48) 


Tatarski  [26]  solved  for  the  inverse  of  Eq.  (2.48)  in  terms  of  the  derivatives  of  the 
structure  function. 


<E'n(«) 


Jo 


°°  sin(«;r)  d 


Kr 


dr 


dr. 


(2.49) 


Equation  (2.49)  can  be  used  to  find  the  power  spectrum  from  a  given  structure 
function.  For  example,  Tatarski  [25]  derived  the  index  structure  function  in  terms 
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of  Kolmogorov’s  inertial  subrange  model 


Dn{r,z)  =  Cl{zy^\ 


(2.50) 


By  substituting  Eq.  (2.50)  into  Eq.  (2.49)  one  can  show 

i,{K,z)=0M3Cl(z)K-''/\ 


(2.51) 


To  find  the  power  spectrum  associated  with  an  arbitrary  structure  function, 


Dn{r,z)  =  Pizy^, 


(2.52) 


one  substitutes  Eq.  (2.52)  into  Eq.  (2.49)  and  evaluates  the  integral.  In  Example 
II. 2  7  was  shown  to  be  7  =  a  —  3,  thus, 


Dn{r,a,z)  =p{zy  I 


(2.53) 


The  power  spectrum  of  the  turbulence  in  terms  of  the  arbitrary  power  law  structure 
function  was  shown  in  Example  II.  1  of  the  previous  section  to  be. 


$„(«,  Z)  =  2“-V(z)(a2  _  5^  g)^-3/2 


a-2 


5— a 


K 


3  <  a  <  5, 


(2.54) 


where  I3{z)  is  the  structure  constant  similar  to  Cl{z)  and  has  units  of  m^““,  r[cc;]  is 
Euler’s  gamma  function,  and  the  power  law  a  is  defined  as  a  =  7  +  3.  The  function 
0(0;)  which  maintains  consistency  between  the  index  structure  function  and  its  power 
spectrum  (See  Eq.  (2.49))  is  found  by  equating  Eq.  (1.5)  and  Eq.  (2.54)  and  solving 
for  a{a): 


a{a)  =  -  5a  +  6)7r-^/2. 


a-2 


5— a 


3  <  a  <  5. 


(2.55) 
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Figure  2.1  The  consistency  function  a{a). 

This  result  is  mathematically  equivalent  to  Eq.  (2.37)  derived  in  Example  II.2  and 
Eq.  (2.57)  below.  The  consistency  function  a{a.)  is  plotted  in  Fig.  (2.1).  Thus,  the 
index  power  spectrum  for  an  arbitrary  power  law  is 

$„(k,  a,  z)  =  a{a)(3{z)K'~°‘.  (2.56) 

Equation  (2.56)  reduces  to  the  Kolmogorov  power  spectrum  (Eq.  (2.51))  if  a  =  11/3. 
Note  that  Beland  [1]  has  derived  an  alternate  but  equivalent  definition  for  a{a) 

T  f-1 

a{a)  =  -(2)°-^-^/^^  3  <  a  <  5.  (2.57) 

^This  problem  was  solved  as  an  example  in  the  previous  section  on  Mellin  transforms. 
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Equations  (2.52)  and  (2.54)  yield  a  very  important  physical  interpretation.  If  the 
inner  and  outer  scale  effects  are  neglected,  the  power  law  a  must  be  greater  than 
3,  because  when  o:  — >  3,  o(a)  0  and  the  turbulence  power  spectrum  vanishes. 
Likewise,  the  index  structure  function  has  no  dependence  on  the  separation  r  (i.e., 
Dn{T,z)  —  I3{z)r'^)^  implying  no  random  turbulence  exists.  We  will  see  in  the  next 
section  that  the  domain  of  a  is  further  constrained  in  the  wave  structure  function 
calculation. 


2.3  The  Wave  Structure  Function 

In  order  to  estimate  the  effects  of  turbulence  on  optical  wave  propagation,  one 
must  solve  the  wave  equation  for  a  random  media, 

V^E  +  -  0,  (2.58) 


where  n  is  a  random  field  of  index  variations  over  the  region  of  propagation.  If  the 
turbulence  is  weak,  homogeneous,  and  isotropic,  the  method  of  small  perturbations 
can  be  used  to  solve  Eq.  (2.58)  [2,  5,  10,  25].  The  method  of  small  perturbations 
yields  the  two  dimensional  spherical  wave  correlation  functions  (in  a  plane  transverse 
to  the  direction  of  propagation)  for  the  log  amplitude  and  phase. 


B^{p)  =  A'K^k'^  j  dz  J  d«;  kJq  sin^ 


k^z{L  —  z) 


2kL 


^^{k,z)  (2.59) 


and 


B^{p)  =  A'K^k'^  dz  dK  kJo 


cos 


k^z{L  —  z) 

2fcl 


#„(k,2),  (2.60) 


where  Jq  is  a  Bessel  function  of  the  first  kind,  order  0,  k  is  the  optical  wavenumber, 
p  is  the  geometrical  separation  between  points  in  the  plane,  and  L  is  the  length  of 
the  optical  path. 
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In  this  section,  the  wave  structure  function  for  non-Kolmogorov  turbulence 
Dy,{p)  is  derived  using  Eq.  (2.59)  and  Eq.  (2.60).  The  wave  structure  function  is 
defined  as  the  sum  of  the  log-amplitude  and  the  phase  structure  functions, 

DM  =  DM  +  DM-  (2-61) 


The  relations 

D^{t)  =  21B^(0)  -  B;,(r)l,  (2.62) 


and 

DM  =  2[5^(0)  -  BMl  (2-63) 

can  be  derived  in  a  fashion  similar  to  Eq.  (2.43).  Using  Eqs.  (2.62)  and  (2.63)  and 
substituting  Eqs.  (2.56),  (2.59)  and  (2.60)  into  Eq.  (2.61)  yields 


D^{p)  =  %'K^k'^a{a) 


K 


1— a 


l-^o(^) 


d«d2:. 


(2.64) 


where  a(a)  is  the  consistency  function,  I3{z)  is  the  index  structure  constant  along 
the  optical  path,  k  is  the  spatial  wavenumber  of  the  turbulent  eddies,  Jq  is  a  Bessel 
function  of  the  first  kind,  order  0,  and  p  is  the  geometric  distance  between  two 
points  in  the  plane  transverse  to  the  direction  of  propagation.  Equation  (2.64)  can 
be  evaluated  using  Mellin  transforms.  Sasiela  [23]  gives  the  necessary  identity, 

M[l  -  Jo(^)]  =  -2''‘ -2  <  »(»)  <  0.  (2-65) 

where  T  is  Euler’s  gamma  function.  Making  the  change  of  variables 


X  = 


Kpz 

~ir 


xL 

K  =  - 

pz 


da?  =  ^dK 
d/c  =  —dx 
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and  letting  s  —  2-  a,  Eq.  (2.64)  becomes 


fiL  /*oo  r  / tTj\^  ^ 

DM  =  sMa{a)  ^  f  —  j  [1  -  Mx)]dxdz.  (2.66) 

Rearranging  terms  yields 

fL  poo 

Dw{p)  =  87r^/!;^a(a)L^““p"“^  /  (3{z)z'^~‘^dz  /  x^~\l  -  Jo(a;)]da;.  (2.67) 

JQ  •'0 

Applying  the  Mellin  transform  in  Eq.  (2.65)  and  remembering  a(a)  is  defined  from 
3  <  CK  <  5, 

2—oc 

DM  =  -(2)^-VA:^a(o)L2-“p“"2r  ^  /  f3{z)z'^-^dz  3  <  a  <  4.  (2.68) 

i/  0 

2  . 

If  P{z)  is  constant  over  the  optical  path,  Eq.  (2.68)  can  be  simplified.  Evaluating 
the  integral  for  this  case  yields 

£)^(p)  = -(2)4-Vfc2a(a)L2-V"^r  ~  3  <  a  <  4.  (2.69) 

'  '  a  a  —  1 

2  . 

Rearranging  Eq.  (2.69)  yields, 

DM  =  -(2)^-“ir^fe^a(a)/; J. j 3  <  a  <  4.  (2.70) 

Equation  (2.70)  is  the  spherical  wave  structure  function  for  turbulence  characterized 
by  an  arbitrary  power  law  a.  To  get  the  plane  wave  structure  function  one  excludes 
the  {a  —  1)  term  in  the  denominator  of  Eq.  (2.70).  Equation  (2.70)  will  be  used 
extensively  in  Chapters  3  and  4  to  compute  the  mutual  coherence  function  and  Strehl 
ratio  of  a  laser  beam. 
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2.4  ^  Generalized  Expression  for  Fried’s  Coherence  Diameter  ro 


Many  of  the  results  in  this  thesis  are  normalized  by  atmospheric  and  path 
conditions.  Fried  [8,  9]  accomplished  the  normalization  by  defining  a  “coherence 
diameter”,  r^,  which  takes  into  account  the  strength  of  the  turbulence  and  the  path 
length.  Fried  [8]  defines  the  plane  wave  structure  function  as 

DM  =  6.88  (2.71) 

Wq  / 

where 

The  arbitrary  constant  6.88  is  defined  more  accurately  in  [8]  as  2(24/5r(6/5))^/®. 
This  constant  was  chosen  to  force  Fried’s  resolution  metric  77^/7?-max  [8]  to  have  a 
certain  asymptotic  behavior.  Since  Fried’s  definition  of  ro  only  applies  to  an  11/3 
power  law  in  the  Kolmogorov  spectrum,  a  more  general  expression  is  presented  here. 

Our  definition  for  the  spherical  wave  structure  function  is  rewritten  from  Eq. 

(2.70) 


DM  =  -(2)^-“7r2fc2a(a)L"-V“-2r 


2-a 

2 

a 

2 


p{z)z'^-^dz  3<a<4.(2.73) 
Jo 


We  need  to  find  a  parameter  po,  which  reduces  to  rg  when  the  power  law  a  11/3, 
such  that 


a-2 


D^{p)  =  Cl  — 

\Po> 


where 


Cl  =  2 


8 


(2.74) 


(2.75) 


—  2  Lq:  —  2j 

The  choice  of  ci  will  become  more  apparent  when  the  asymptotic  behavior  of  Fried’s 
'R-l'Ernax  is  discussed  in  Chapter  4.  Setting  Eq.  (2.73)  equal  to  Eq.  (2.74)  and 
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solving  for  po  yields 


Po 


Ci{“  -  I)r  [f 

-(2)4-“7r2fc2a(a)L2-«r 

p{z)z<^-^dz 

1 

ex-2 


3  <  a  <  4.  (2.76) 


If  p(z)  is  constant  over  the  path  then  Eq.  (2.76)  reduces  to 


Po 


[  Ci(a  -  I)r  [; 

^1  1 

1  — (2)^““7r2fc2a(a)r 

'2-a' 

2 

pL] 

3  <  cc  <  4. 


(2.77) 


Again,  for  the  plane  wave  expression  of  po  the  {a  —  1)  term  is  omitted. 

Figure  (2.2)  is  a  plot  of  po  as  a  function  of  a  for  various  values  of  /?L.  It 
is  important  to  note  that  po  scales  with  In  Section  2.2,  it  was  shown  that 


pz,  =  4  X  lO”'* 
PL  =  4x10'“' 
PL  =  4x10"” 


Figure  2.2  The  generalized  coherence  diameter  po  as  a  function  of  a  normalized  by 
PL.  The  wavelength  is  1.3  microns. 
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turbulence  with  power  laws  o:  <  3  can  not  exist,  because  the  index  structure  function 
has  no  dependence  on  r  when  a  =  3.  This  is  reflected  in  Fig.  (2.2)  as  well,  since 
Po  oo  when  a  — >  3.  When  o:  — ^  4  we  see  >  0,  this  can  be  explained  by 
the  square  law  dependence  of  the  wave  structure  function  (Eq.  (2.74)  when  a  =  4. 
Heidbreder  [12]  showed  that  the  phase  structure  function  in  the  near  field  limit,  which 
is  equivalent  to  the  wave  structure  function  presented  here,  induces  pure  random  tilt 
on  the  wavefront.  Since  a  pure  random  tilt  shifts  the  beam  olf  axis  for  each  realization 
of  turbulence  (short  exposure)  the  average  irradiance  (long  exposure)  will  be  the 
average  of  the  unabberated  beam  shifted  all  over  the  observation  plane.  This  results 
in  a  wide  spot  with  a  very  low  irradiance  peak.  This  is  exactly  what  is  happening 
in  Fig.  (2.2).  The  worst  turbulence  to  have  in  an  uncompensated  beam  projection 
system  is  turbulence  whose  structure  function  only  induces  tilt  on  the  wavefront, 
thereby  shifting  its  short  exposure  spot  continually  around  the  observation  plane. 
This  eflPect  will  be  discussed  further  in  Chapter  4. 

2.5  The  Variance  of  the  Log  Amplitude  Fluctuations  and  Limitations  of  the  Theory 

Using  the  method  of  small  perturbations  to  solve  the  wave  equation  (See  Eq. 
(2.58))  is  limited  to  regimes  where  the  tubulence  is  weak  and  the  propagation  dis¬ 
tance  is  short.  For  Kolmogorov  turbulence  the  method  of  small  perturbations  has 
been  shown  to  be  valid  when  the  log  amplitude  variance  is  small  (for  cr^  <  0.3  -  0.5) 
[5].  When  the  log  amplitude  variance  is  greater  than  0.3  —  0.5  a  multiple  scatter¬ 
ing  theory  is  required.  For  purposes  of  this  discussion,  it  is  assumed  the  criterion 
cr^  <  0.3  —  0.5  holds  for  propagation  in  an  arbitrary  power  law  turbulence  spectrum. 
To  prove  the  relation  cr^  <  0.3  -  0.5  holds  for  an  arbitrary  power  law  turbulence 
spectrum  requires  comparison  with  a  multiple  scattering  theory  which  is  beyond  the 
scope  of  this  thesis. 

One  can  examine  the  variance  of  the  log  amplitude  perturbations  in  order 
to  understand  the  bounds  the  method  of  small  perturbations  places  on  turbulence 
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strength  and  range  as  the  power  law  of  the  turbulence  spectrum  varies.  Here  the 
bounds  are  found  by  constraining  the  log  amplitude  variance  to  0.1  -  0.5  and  solv¬ 
ing  for  the  system  and  atmospheric  parameters  of  interest.  From  the  plane  wave 
correlation  function  [2], 


fL  f°°  ^  .2 

k^{L  —  z) 

dz  dK  K  sin 

Jo  Jo 

2k 

(2.78) 


Substituting  the  turbulence  spectrum  given  by  Eq.  (2.56)  into  Eq.  (2.78)  yields 

k?{L  —  z) 


fL  /*oo 

al  =  f  dzl3{z)  d/c/c^""sin^ 

^  Jo  Jo 


2k 


(2.79) 


Letting  s  =  2  —  a  and  w  =  gives 


7^  =  4:'!r'^k'^a{a)  j  dz  P{z)  dKK^  ^  sin^ 


WK 


(2.80) 


Making  the  change  of  variable  x  =  ^/wK  and  da;  =  ^/wdn  yields 


<y‘ 


=  47r^fc^a(Q:)  j  dz  P{z)  J 


°°  da;  /  X  ^ 
y/w  \y/w ^ 


3  —  1 


sin^  .  (2.81) 


Simplifying  Eq.  (2.81)  gives 


=  47r^fc^a(a)  j  dzw  */^/?(z)  ^  dx  (x)*  ^  sin^ 


X" 


(2.82) 


Equation  (2.82)  is  the  Mellin  transform  of  sin^(x2).  The  applicable  transform  pair  is 


M[sm^{x'^)]  -^r 


s/4 

1/2 -s/4 


-4  <  K(s)  <  0, 


(2.83) 
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thus, 


a^-  -4Tr‘^k'^a{a)  dzw 


2-a 

4 

a 

4 


2  <  0!  <  6. 


Substituting  for  w, 


0(.)V 


2-a 

4 

a 

4 


2  <  a  <  6. 


Simplifying  Eq.  (2.85)  and  remembering  a{a)  is  defined  on  3  <  o;  <  5  gives 


^2  ^ 


2-a 

4 

a 

4 


/  dz^iz)  {L  —  z)  3  <  a  <  5. 
Jo 


If  p{z)  is  constant  along  the  path,  then  we  have 


2  2^ 'K^l'^k  2°‘a{a). 


a 


2-a 

4 

a 

4 


fJL^  3  <  a  <  4. 


In  order  to  generalize  for  various  system  and  atmospheric  conditions, 
multiply  Eq.  (2.87)  by  Po~'^IPo~^  and  rearrange  the  terms  to  yield 


0  —  6  1  /O  -1  2— or  _  01—2 

2  2  2  ci'K^'^k  2  L  2 

X  ru  nO. — 2 


2-a 

4 

a 

4 


a 

2 

2-a 

2 


Define  the  Fresnel  number 

FN  =  VTl, 


(2.84) 


(2.85) 


(2.86) 


(2.87) 

we  can 


(2.88) 


(2.89) 
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where  A  is  the  wavelength  of  the  beam  and  L  is  the  length  of  the  optical  path. 
Substituting  the  definition  of  Fresnel  number  into  Eq.  (2.88)  yields 


a 


2 

X 


3- 

CiTT  2 

4q! 


(2.90) 


Figure  (2.3)  is  a  plot  of  the  log  amplitude  variance  as  a  function  of  power  law 
normalized  by  (FiV/po)""^.  Care  must  be  used  in  interpreting  this  plot.  Remember 


Figure  2.3  The  log  amplitude  variance  as  a  function  of  a  for  normalized  by 
{FNlpor-^. 

for  constant  turbulence  and  path  conditions  po  varies  considerably  with  a  (See  Fig 
(2.2)).  A  key  observation  can  be  made  at  this  point.  When  the  spectrum  contains 
fewer  eddies  of  high  wavenumber,  corresponding  to  higher  power  laws  (a  — >  4), 
scintillation  effects  are  reduced.  This  means  phase  effects  dominate  as  a:  ^  4, 
supporting  the  conclusions  of  Heidbreder.  Figure  (2.4)  contains  plots  of  the  minimum 
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Figure  2.4  The  minimum  po  normalized  by  Fresnel  number  FN  for  the  conditions 
where  =  0.1, 0.3, 0.5. 

{po/FN)°‘~‘^  which  can  be  examined  using  this  theory  for  various  values  of  cr^.  One 
can  reject  errors  due  to  multiple  scattering  effects  by  restricting  use  of  the  theory 
presented  in  this  thesis  to  the  areas  above  the  curves.  The  area  below  the  curves  is 
considered  “strong”  turbulence  and  the  theory  presented  here  no  longer  applies. 
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III.  The  Extended  Huygens- Fresnel  Principle 

In  this  chapter,  the  extended  Huygens- Fresnel  principle  is  used  to  derive  the 
mutual  coherence  function  for  an  arbitrary  wave  in  a  vacuum  and  in  turbulence. 
Once  the  MCF  for  an  arbitrary  wave  in  turbulence  is  found,  it  is  a  simple  matter  to 
determine  the  intensity  at  any  point  in  the  observation  plane.  The  on-axis  intensity  is 
of  particular  interest,  because  it  allows  one  to  easily  generate  a  useful  figure-of-merit; 
the  Strehl  ratio  {STZ).  The  Strehl  ratio  can  be  shown  to  be  the  ratio  of  the  on-axis 
intensity  of  a  laser  beam  in  turbulence  to  that  of  the  laser  beam  in  a  vacuum.  Thus, 
the  Strehl  ratio  provides  a  useful  metric  for  examining  the  impacts  of  turbulence  on 
laser  beam  propagation.  In  this  chapter,  the  general  equations  for  the  MCF  and 
intensities  are  derived.  In  Chapter  4,  the  equations  for  the  on-axis  intensity  of  a 
general  wave  are  used  to  calculate  the  Strehl  ratio  of  a  focussed  constant  amplitude 
beam. 

3. 1  The  Mutual  Coherence  Function  for  Quasimonochromatic  Waves  in  a  Vacuum 

Following  the  approach  of  Lutomirski  and  Buser  [17],  the  MCF  for  quasi¬ 
monochromatic  light  is  defined  as 

MCF(Pi,P2)  =  {U{xi,yi)U*{x2,y2)) ,  (3.1) 

where  pi  and  p2  are  points  in  the  observation  plane,  (xi,yi)  and  {x2,y2)  are  the 
X  and  y  coordinates  of  pi  and  p2  and  the  angle  brackets  denote  the  statistical 
expectation  operator.  Figure  (3.1)  shows  the  geometry  of  the  propagation  scenario. 
From  the  Huygens- Fresnel  principle 

U{x,y)  =  jjG{x,y,^,r})UA{^,v)d(dr],  (3.2) 


3-1 


where  ^  and  ry  are  the  coordinates  of  a  point  in  the  source  plane  and  the  Green’s 
function  for  propagation  of  a  spherical  wave  is 


Ndw  wavefront 


Turbulent  Atmosphere 


Image  Plane 


Figure  3.1  A  wave  propagating  from  the  {C,v)  plane  to  the  {x,y)  plane,  using  the 
Huygens- Fresnel  principle.  The  primary  wavefront  is  the  summation  of 
the  envelopes  of  spherical  wavelets  whose  centers  where  located  on  an 
earlier  wavefront. 


G{x,y,^,r])  = 


exp  [jk{{x  -  +  (y  -  vf  + 


((a;  -  0^  +  {y-  nf  + 


(3.3) 


The  integrals  in  Eq.  (3.3)  are  evaluated  over  the  source  aperture.  Using  the  paraxial 
approximation  {z^  »  {x-i)'^-{-{y-rif),  the  numerator  of  Eq.  (3.3)  is  approximated 

by 


exp  \jk  ((a;  -if  ^ 

and  the  denominator  by 


1/2- 


exp 


\  2Z  ; 


(3.4) 


((i  -if  +  (v-r]f  +  -*  z 


(3.5) 


Substituting  Eqs.  (3.4)  and  (3.5)  into  Eq.  (3.2)  yields, 


,  ,  —jk 

U{x,y)  =  —  exT 


2'n'z 


ff 

L77(^-0'  +  (y-^)'Y 

r  V  2-  )\ 

UA{iyV)<ii^V-  (3-6) 
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Expanding  the  terms  of  the  exponential  and  rearranging,  the  field  is 
U{x,y)  = 

^exp[ifcz]exp  + 

X  yy  exp  +  exp  -^{xi-^yrj)  Ua{^,v)^C^V-  (3-7) 

Eq.  (3.7)  is  the  standard  Fresnel  diffraction  formula.  Substituting  Eq.  (3.7)  into  the 
definition  of  the  MCE,  Eq.  (3.1),  yields 

MCF{xi,yi,X2,y2,z)  = 

X  yy  exp  Yz^il  +  nl)  exp  ^^(xi^i  +  yi77i)  UA{^i,Vi)d(idyi 
X  yyexp  -^{^l  +  vD  exp  ^(^26  +  2/2^2)  UX{^2,V2)d^2dy2^  .  (3.8) 

Since  the  wave  is  propagating  in  a  deterministic  fashion,  the  expectation  operator 
can  be  neglected  and  the  MCE  for  a  coherent  wave  in  a  vacuum  is 

MCF{xi,yi;x2,y2,z)  = 

^  //// 

r _ 

X  exp  — ^(xi^i  +  yiyi  +  ^26  +  y2’?2) 

X  Ua{Ci,Vi)UX{C2,  ?72)d6dmd6d/72-  (3.9) 

By  defining 

ijj  =  exp  ^(x^  +  yl)  -  (x^  +  yl)  ,  (3.10) 
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the  MCF  for  a  quasimonochromatic  wave  propagating  in  a  vacuum  is 


MCF{xi,yuX2,y2yz)  = 

r _ 7*Ai 

X  exp  —^{xi^i  +  yiTii  +  X2^2  +  2/2^72) 

X  Ua{^i,Vi)UXU2,  ^?2)d6dmd6d772-  (3-11) 


3.1.1  The  irradiance  of  a  quasimonochromatic  wave  in  a  vacuum.  The 
irradiance  at  a  point  {xi,y\)  can  be  determined  from  the  MCF  by  letting  X2  Xi 
and  y2  y\  in  Eq.  (3.11), 


I{xi,yi,z)  =  MCF{xi,yx\Xx,yi,z)  = 


jk 

2z 


m  +  vD-i^l  +  vD) 


X  exp 


-jk 


(a;iCi  +  yiVi  +  +  yi’72) 


X  UAi^uVi)Ul{C2,  ?72)d6d?7id6d?72- 


(3.12) 


If  we  consider  only  the  on-axis  irradiance  xi  — >  0  and  yx  — ^  0,  Eq.  (3.12)  becomes 


J(0,0,z)  = 

1=17/// 


exp 


^^{{Cl  +  vl)-{f2+vl)) 


X  Ua{^i,Vi)Ua{^2,  '?72)d6dr?id6d772- 


Equation  (3.13)  is  a  separable  integral, 


(3.13) 
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J(0,0,z)  = 

k  .o  f  r  \jk 


// 


exp 


zitif2 

2z 


t^A(6,»7i)d6d?7i 

J 

i^l  +  vl)  C^I(6,%)d6d?72. 


(3.14) 


5.^  The  Mutual  Coherence  Function  for  Quasimonochromatic  Waves  in  Non- Kolmogorov 
Turbulence 

In  this  section  the  MCF  for  an  arbitrary  wave  in  non- Kolmogorov  turbulence 
is  derived.  The  derivation  proceeds  as  in  the  last  section,  except  for  the  substitution 
of  a  different  Green’s  function. 

The  MCF  for  quasimonochromatic  light  is  defined  as 


MCF{pi,P2)  =  {U{xi,yi)U*{x2,y2)),  (3.15) 


where  pi  and  p2  are  points  in  the  observation  plane,  and  (xi,yi)  and  {x2,y2)  are 
the  X  and  y  coordinates  of  pi  and  p2. 

From  the  Huygens-Fresnel  principle; 

U{x,y)^-^  JJ  G{x,y,f,T])UA{^,v)dCdr],  (3.16) 


where  ^  and  rj  are  the  coordinates  of  a  point  in  the  source  plane  and  G{x,  y,  y)  is 
the  Green’s  function  for  propagation  in  turbulence.  The  integrals  are  again  taken 
over  the  source  aperture.  If  we  define  G'{x,y,(,ri)  as  the  ratio  of  G{x,y,(,r])  to  its 
value  in  the  absence  of  turbulence  [17], 


exp[jfc((a;  -  ^)^  +  (y  -  vf  + 

((a:  -  ^)2  +  (y  -  77)2  -1-  z2)i/2 


G\x,y,i,y). 


G{x,y,^,y)  = 


(3.17) 


Applying  the  paraxial  approximation  and  substituting  Eq.  (3.17)  into  Eq.  (3.16) 
yields 


U{x,y) 


-jk 

2‘kz 


exp[jkz] 


X  jf  G'(i>  »,{.>?)  exp 


-  0"  +  (v  -  ’>)' 


UA{i,ri)d(d,,.  (3.18) 


Substituting  Eq.  (3.18)  into  Eq.  (3.15), 


MCF{xi,yi]X2,y2,z)  = 

Jjj f  G'{^uyu^i^Vi)G'*{x2,y2,$2,m) 

^  ^ {xi  -  +  {yi  -  mf' 


X  exp  \  jk 


2z 


X  exp 


-jk 


'  ix2  -  ^2f  ±(jj2-V^ 


2z 


X  UA{^um)UX{^2,V2)d$idriid$2<iV2 


(3.19) 


Interchanging  the  order  of  integration  and  expectation  and  defining 


H{xi,yi,^i,r]i,X2,y2,(2,V2,z)  =  {G' {xi,yi,  z)G'* {x2,y2,  C2,'n2,  z)),  (3.20) 


yields 


MCF{xi,yi,X2,y2,z)  = 

J fJJ  H{xi,yi,^i,yi,X2,y2,C2,V2,z) 

X  exp  j^((^i  “  +  (^1  “  ^i)^)  ~  ((®2  -  6)^  +  (y2  -  V2f)) 

X  Ua{^i,Vi)UX{(2,  m)d^idmd^2dri2.  (3.21) 
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Expanding  the  terms  in  the  exponential  and  rearranging,  Eq.(3.21)  becomes 


MCF{xuyi]X2,y2,z)  = 

^  Jill 

X  exp  [  +  vl)  -  (^2  +  vl)  +  2a:26  -  2xi^i  +  ^y2V2  -  ‘^yiVi) 

\_2z 

X  UaUi,  Vi)U*a{$2,  ^2)d6dr?id6d772. 


(3.22) 


By  defining 


V'  =  exp  ^(cci  +  yl)  -  {xl  +  yj)  , 


(3.23) 


the  MCE  of  a  quasimonochromatic  wave  in  non-Kolmogorov  turbulence  is 


MCF{xi,yi]X2,y2,z)  = 

(^)V  ////  H{xi,yi,^uyi,X2,y2,(2,V2,z) 

X  exp  |^^((^i  +  ^i)  -  (Ci  +  vl)  +  2a;26  -  2a;i^i  +  2y2y2  -  2yiyi) 
^Ua{^i,Vi)UX{^2,  ?72)d6d77id^2d772-  (3-24) 

Equation  (3.24)  is  the  major  result  of  this  chapter.  We  will  see  in  the  next  chapter 
that  H{xi,  yi,  6,  ?7i,  ®2,  ^2, 6,  V2,  z)  takes  the  form 

H{xi,yi,^uyi,X2,y2,^2,V2,z)  =  exp  -^DM  (3-25) 

where  Dyj{p)  is  the  wave  structure  function  derived  in  Chapter  2  and  p  is  the  geo¬ 
metric  separation  of  the  points  (Ci,?7i)  and  (^21^2)  iii  fh®  source  plane. 

3.2.1  The  irradiance  of  a  quasimonochromatic  wave  in  non-Kolmogorov  tur¬ 
bulence.  The  irradiance  at  a  point  (x,  y)  in  the  observation  plane  can  be  deter- 
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mined  by  letting  X2  — *■  Xi  and  ?/2  — ^  Vi  iii  Eq.  (3.24), 


(— )^y' IJj  H{xi,yi,^i,rji,X2,y2,^2,V2,z) 

X  exp  +  yl)  -  +  vl)  +  2x1^2  -  2xi^i  +  2yiy2  -  2yiyi) 

2z 

xUa{^i,Vi)UaU2,  %)d6d^id6d?72- 


(3.26) 


Considering  only  the  on-axis  irradiance,  let  rci  — ^  0  and  yi  — >  0,  then 


7(0,0,  .2)  = 

I jjj  H{xi,yx,iuVuX2,y2,i2,m,z) 


X  exp 


j_k 

2z 


X  UaHi  ,  Vi)Ul{^2,  V2)d^idvid^2dv2- 


(3.27) 


We  will  use  the  Eqs.  (3.14)  and  (3.27)  for  the  on-axis  intensities  in  vacuum  and 
in  turbulence  extensively  in  Chapter  4  to  compute  the  Strehl  ratio  of  a  focussed, 
constant  amplitude  beam. 
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IV.  The  Strehl  Ratio  of  a  Laser  Beam  in  Non- Kolmogorov 
Atmospheric  Turbulence 

In  order  to  calculate  the  impact  turbulence  has  on  laser  beam  propagation, 
one  can  use  many  different  performance  metrics.  A  common  performance  metric  is 
the  Strehl  ratio.  The  Strehl  ratio  is  defined  as 


sn  = 


/(0,0) 


in  turbulence 


7(0,0) 


in  vacuum 


(4.1) 


where  7(0,  0)  is  the  on- axis  irradiance. 

In  this  chapter  the  on-axis  irradiances  of  a  wave  in  a  vacuum  and  in  turbulence, 
hence  the  Strehl  ratio,  are  calculated  using  the  mutual  coherence  function  (MCF) 
derived  in  Chapter  3  and  the  wave  structure  function  D^ip)  derived  in  Chapter  2. 
In  Section  4.2,  the  Strehl  ratio  as  a  function  of  power  law  a  is  plotted  for  constant 
values  of  the  generalized  coherence  diameter  po-  Next,  the  asymptotic  behavior  of 
Fried’s  resolution  metric  7^/7^^lax  is  examined  [8].  It  is  from  this  behavior  that  the 
constant  Ci  in  Eq.  (2.74)  is  determined.  Additionally,  T^/TT-max  is  shown  to  be  the 
Strehl  ratio  multiplied  by  [Dj poY  or 


n 

TT^max 


^sn 


(4.2) 


Equation  (4.2)  is  plotted  for  various  power  laws  in  Section  4.3. 

Great  care  must  be  used  in  finding  a  reasonable  basis  on  which  to  compare  the 
Strehl  ratios  of  beams  propagating  through  different  power  spectra.  Consider  the 
two  power  spectra 


$,,,(«, ai,^o)  =  a(ai)^o«  a2,/5o)  =  a(o!2)/5oK 


(4.3) 
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If  cti  does  not  equal  0:2,  then  the  power  in  any  finite  bandwidth  of  the  two  spectra 
will  not  be  equal.  Thus,  in  a  plot  of  the  Strehl  ratio  as  a  function  of  the  power  law 
a,  it  is  difficult  to  separate  the  effects  of  the  total  power  in  the  turbulence  from  how 
that  power  is  weighted  in  spatial  wavenumber.  We  shall  see  in  this  chapter  that, 
the  choice  of  basis  for  comparing  different  power  laws  greatly  impacts  the  trends  of 
Strehl  ratio. 

The  latter  sections  of  this  chapter  contain  plots  of  the  Strehl  ratio  as  a  function 
of  the  power  law  under  different  bases  of  comparison.  Ultimately,  any  individual 
basis  boils  down  to  a  rule  on  how  (3L  is  chosen  as  the  power  law  varies.  First, 
PL  is  fixed  and  the  power  law  a  is  allowed  to  vary  over  the  domain,  3  <  a  <  4. 
Following  this  comparison,  PL  is  chosen  such  that  the  power  spectra  are  constant  at 
a  specific  wavenumber  Kq-  Next,  pL  is  chosen  to  be  such  that,  the  power  over  a  finite 
bandwidth  is  constant  as  the  power  law  varies.  We  will  see  that  the  results  of  the  two 
previous  comparisons  are  similar.  Next,  PL  is  chosen  such  that,  the  piston-removed 
wave  variance  over  the  receiver  aperture  is  constant  as  the  power  law  varies. 


4.1  Strehl  Ratio  for  a  Focused  Constant-Amplitude  Beam 

In  this  section  an  equation  for  the  Strehl  ratio  is  derived  using  the  wave  struc¬ 
ture  function  derived  in  Chapter  2  and  the  extended  Huygens-Fresnel  principle  de¬ 
scribed  in  Chapter  3.  First,  let  the  field  in  the  source  aperture  be  a  constant- 
amplitude  beam  focussed  at  distance  z. 


U{^,T])  =  exp 


j_k 

2z 


(4.4) 


4.1.1  Evaluation  of  the  on-axis  irradiance  in  turbulence.  Substituting  Eq. 
(4.4)  into  Eq.  (3.27)  yields 
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7(0, 0,z) 
k 


in  turbulence 


J f f f  H{xuyu(i,VuX2,y2, ^2, V2,  z)d^idr]idC2dV2-  (4.5) 

Lutomirski  and  Buser  [17]  show  that  H{xi,yi,^i,r]i,X2,y2,^2,y2,  z)  in  Eqs.  (3.24), 
(3.26)  and  (3.27)  is 


2/i)  ^2)  2/2)  ^2)  ^2)  .2)  exp 


-  6)^  +  {vi  -  V2Y)^^^  )  (4.6) 


where  is  the  spherical  wave  structure  function  derived  in  Chapter  2  and  p  — 
((^1  —  6)^  +  (^1  ■"  '^2)^)^^^  is  the  geometric  distance  between  two  points  in  the  source 
plane.  Substituting  Eq.  (4.6)  into  Eq.  (4.5)  yields 


7(0,  0,  2^)in  turbulence 

'If  III! 


exp 


■DM 


d^idr]id^2<iV2- 


(4.7) 


Equation  (4.7)  is  a  four  dimensional  integral  over  the  source  aperture.  We  define  the 
circular  aperture  function 

1  \^d-^Viv\  <  D/2 

0  otherwise, 


where  |  and  fj  are  unit  vectors  in  the  source  plane  and  D  is  the  aperture  diameter. 
Substituting  Eq.  (4.8)  into  Eq.  (4.7)  yields 


5  ^Jin  turbulence 

jU  /*oo  roo  1*00  roo  ^  y. 

/  /  /  w(u  +  nifi)w(U  +  rnfi) 

Z'KZ  J—ooJ—ooJ—ooJ—00 

X  exp  d6dmd6d?72.  (4.9) 
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Equation  (4.9)  is  converted  to  vector  notation  with  the  following  change  of  variables 


thus, 


r  =  (6  -  6)^^  +  {Vi  -  V2)fj  dr  =  d^idr^i 
^  +  6)1  +  (^1  +  V2)v)  dr'  =  d6d772, 

f  +  |(f)  =  +  r}i-n  F  -  1(f)  =  +  V2V- 


Substituting  Eqs.  (4.10)  and  (4.11)  into  Eq.  (4.9)  yields 


(4.10) 


(4.11) 


Oj  turbulence 

U  /*oo  roo  1 

X  exp  [--E>-«,(|fl)  dfaf. 


5(?)) 


(4.12) 


Because  hE(f'  +  |(f))  and  l^(f  -  |(f))  represent  circular  apertures  they  are  inde¬ 
pendent  of  orientation  in  the  source  plane,  thus,  Eq.  (4.12)  can  be  written  as 


f(0,  0,  turbulence 

rji  ”'<1^  -  jwi) 

X  exp  dfdf'. 


dfdr . 


(4.13) 


The  integral  over  f'  can  be  evaluated  in  closed  form  [9], 


I  W{\f  +  icr^DH/d?  -  i(f)|)dr  = 


(4.14) 


otherwise. 
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where  D  is  the  diameter  of  the  aperture.  Substituting  Eq.  (4.14)  into  Eq.  (4.13) 
yields 


I (Oj  0)  ^)in  turbulence 

X  exp  j--i:)^(lf|)  dr. 

Changing  to  polar  coordinates,  the  on-axis  irradiance  becomes, 


(4.16) 


/(0,0,2)i„r„,b„,ence  =  (^)^27rj^  TsfC^r) cxp  -^D^{r)  dr,  (4.16) 


where, 


Koir)  =  i  D^cos  l(^)  -  r{D^  - 


(4.17) 


4.1.2  Evaluation  of  the  on-axis  irradiance  in  a  vacuum.  The  on-axis  irra¬ 
diance  in  vacuum  is  given  by  Eq.  (3.14) 


/(0,0,4 

r  k 


+  U*j^{C2,V2)d^2dV2- 


(4.18) 


Substituting  the  field  (Eq.  (4.4))  and  the  aperture  functions  (Eq.  (4.8))  and  evalu¬ 
ating  the  integrals  yields 


7(0,0,r)..._  =  (ffj 


(4.19) 
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Using  the  definition  of  Strehl  ratio  we  obtain 


57^  = 


/(0,0)i„ 


turbulence 


/(0,0), 


in  vacuum 


dr. 


(4.20) 


Equation  (4.20)  is  evaluated  numerically  in  the  latter  sections  of  this  chapter,  as 
the  wave  structure  function  is  varied  according  to  a  rule  on  how  f3L  is  chosen  for  a 
particular  power  law  a. 


4.2  Asymptotic  Behavior  of  Fried’s  Resolution  Metric  for  Non- Kolmogorov  Turbu¬ 
lence 

In  the  previous  section,  the  Strehl  ratio  was  found  to  be 


<57e  = 


7(0,0), 


in  turbulence 


7(0,0), 


in  vacuum 


dr. 


(4.21) 


Equation  (4.21)  can  be  simplified.  Making  the  change  of  variables 


r 

u  =  — 
D 

dr 


d«  = 


D 


r  =  Du 
dr  =  Ddu, 


(4.22) 


and  substituting  Eq.  (4.17)  into  Eq.  (4.21)  yields 


572.  =  -^  [  -(D^cos  ^{u)  -  Duy/D"^  -  D'^u‘^)exp  -l-Dy,{Du)  D^udu.  (4.23) 
Jo  2  L  2  . 


Substituting  Eq.  (2.74)  into  Eq.  (4.23)  and  simplifying  yields 


STZ  =  —  /  -(cos  ^(u)  —  uVl  —  u^)exp 
TT  Jo  2 


1 

(D 

-0^1 

— u 

2 

\Po  ) 

udu. 


(4.24) 


In  order  to  evaluate  Eq.  (4.24)  we  need  to  choose  an  appropriate  constant  ci.  This 
will  be  done  in  the  same  manner  as  Fried  [8].  Fried  chose  the  constant  ci  based  on 
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the  asymptotic  behavior  of  a  performance  metric  similiar  to  the  Strehl  ratio.  He 
called  this  metric  the  resolution  IZ.  The  resolution  TZ  is  defined  as 


7^  =  27r 


(4.25) 


where  H{oj)  is  the  long  exposure  optical  transfer  function  (OTF)  of  an  incoherent 
imaging  system  in  the  presence  of  turbulence.  Fried  examined  the  ratio  TZlTZmaxi 
where  TZmax  is  defined  as  the  resolution  of  an  imaging  system  in  turbulence  with  an 
infinite  aperture  diameter 

Umax  =  lim  n.  (4.26) 

D— voo 

The  ratio  of  'IZlTZ^ax  yields  an  expression  similar  to  Strehl  ratio  in  Eq.  (4.24). 


TZ/n^ax  =  /  J(cos  ^(u)  -wVl-'u2)exp  -ici 

irri  JQ  2  2  kTo  / 


1  fD 

u 
r. 


udu,  (4.27) 


where  r^,  is  the  atmospheric  coherence  diameter  for  Kolmogorov  turbulence  (See  Eq. 
(2.72))  [8].  Thus,  Fried’s  resolution  metric  UjUmax  for  an  arbitrary  power  law  can 
be  found  by  multiplying  Eq.  (4.24)  by  {D/poY 


TZlUmax  =  Sn{Dlpof  = 


32£>- 


■  j ^l{cos-^{u)  -  U\/1  -  «2)exp 


udu. 


(4.28) 


To  find  the  constant  Ci,  make  the  change  of  variables 


x  =  ^  u  =  ^ 

Po  D 

dx=^  du  =  ^ 

Po  ^ 


Thus,  Eq.  (4.28)  becomes 
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toil-* 


'R'/T^max  — 

32  /-f-  1  /  _i  (PoX\  (poX\  f  fpoXV 

7/0  2^  {-D)-[-D)r-[-D) 

r  1  -1 

£cda;. 


X  exp 


a-2 


Let  the  aperture  diameter  increase  until  D 


(4.29) 


00,  thus 


{cos-^{poxlD)  -  {poxlD)^Jl  -  {pox/Df^j  7r/4  and  Eq.  (4.28)  takes  the  form 


foa  pi 

lim  {n/Tlmax}  =  8  /  exp  --Cia;“" 
D->oo  Jo  L  2 


xdx. 


(4.30) 


By  making  another  change  of  variables 


*^1  a-2 

V  =  —X 

2 


X 


=  (!)"■’ 

dv  =  (jj.  — 


(4.31) 


Eq.  (4.30)  reduces  to 

2  ^ 

lim  {Rln^ax}  = - 7  f— '  /  expi-riJu^-Mu. 

D^oo  a  — 2  Vci/  ^0 


(4.32) 


The  integral  in  Eq.  (4.32)  can  be  recognized  as  r[2/a;  -  2],  thus  Eq.  (4.32)  is 

9 

2 


8  /  2  \  ^ 

lim  {RlTlmax]  = - 7  (  — )  r 

z>-»oo"  a  — 2\ci/  Lex  —  2J 


Since  lim^j^oo  {^/^max}  =  TlmaxI'Rmax  =  1  solving  for  Cl  yields 


(4.33) 


Cl  =  2 


a-2  Lcx-2J 


(4.34) 
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If  we  let  a  =  11/3,  Ci  reduces  to  the  constant  in  Fried’s  structure  function 


(4.35) 


To  show  the  behavior  of  an  imaging  system  when  the  aperture  size  is  much  smaller 
than  po  let  D/po  ^  0,  the  exponential  term  in  Eq.  (4.28)  goes  to  1  and 


‘30/12  1  _ 

Tll'R-max  = - ^  /  :^(cos~^(u)  -  uy/l  -  u^)udu. 

'Kpl  Jo  2 


(4.36) 


Equation  (4.36)  reduces  to 


16  /  D\^  fTT  TT  \ 

— (g  le)- 


-) 

PoJ 


(4.37) 


Figure  (4.1)  shows  TllTlmax  as  a  function  of  Djpo  for  various  power  laws. 
Notice  the  only  effect  the  power  law  has  on  TllTlmax  is  to  round  off  the  knee  of  the 
curve.  This  is  another  result  of  the  definition  of  po-  Remember  po  varies  greatly  with 
the  power  law  for  a  constant  integrated  turbulence  strength. 

4.3  Comparing  Strehl  Ratios  in  Arbitrary  Power  Law  Turbulence  Spectra  when 
D I  Po  is  Constant 

Figure  (4.2)  is  a  plot  of  the  Strehl  ratio  as  a  function  of  the  power  law  a  for 
various  values  of  D  j po-  The  Strehl  ratio  is  computed  using  Eq.  (4.24).  Notice  that 
the  Strehl  ratio  does  not  seem  to  change  appreciably  with  power  law.  Remember 
that  for  a  fixed  integrated  turbulence  strength  jTL,  po  varies  greatly  with  the  power 
law  a  (See  Fig.  (2.2)).  Thus,  in  Figure  (4.2)  D  is  varying  just  as  greatly  with  the 
power  law  to  maintain  a  constant  ratio  of  DJ po-  It  is  clear  that  another  point  of 
view  must  be  taken  to  understand  how  the  Strehl  ratio  varies  with  power  law. 
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Figure  4.1  TllUmax  as  a  function  of  Dfpo  for  various  power  laws  (See  Eq.  (4.28)). 

4.4  Comparing  Strehl  Ratios  in  Arbitrary  Power  Law  Turbulence  Spectra  when  (3L 

is  Constant 

The  next  basis  on  which  to  compare  Strehl  ratios  as  the  power  law  a  varies  is  to 
leave  I3L  constant.  This  basis  is  appropriate  for  the  ARGUS  anemometer  case  where 
P  may  be  estimated  correctly,  but  the  uncertainty  in  the  actual  value  of  the  power 
law  is  large  due  to  the  unknown  anemometer  probe  frequency  response.  Figure  (4.3) 
is  a  plot  of  the  Strehl  ratio  as  a  function  of  the  power  law  a  for  a  1  m  diameter, 
circular  source  aperture.  The  wavelength  is  1.3  microns.  Notice  Fig.  (4.3)  follows 
the  same  general  trends  as  po  in  Fig.  (2.2).  We  would  expect  this  behavior  since  f3L 
is  fixed  in  both  cases.  The  Strehl  ratio  approaches  1  when  a  =  3,  because  as  was 
shown  in  Chapter  2,  the  turbulence  vanishes.  The  Strehl  ratio  approaches  0  when 
a  =  4  because  the  abberations  induced  on  the  optical  wave  are  pure  random  tilts  in 
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Figure  4.2  The  Strehl  ratio  as  a  function  of  a  {See  Eq.  (4.24).  Dj po  is  fixed  at  the 
value  specified  in  the  legend. 

the  source  plane.  These  random  tilts  cause  a  greatly  increased  long  exposure  spot 
size  due  to  beam  wander.  Note:  this  interpretation  is  only  valid  if  inner  and  outer 
scale  effects  are  neglected. 

4.5  Comparing  Strehl  Ratios  in  Arbitrary  Power  Law  Turbulence  Spectra  when  (3L 
is  such  that  a{a)f3LK~°'  is  Constant 

In  this  comparison  fJL  is  chosen  such  that  the  power  spectrum  for  each  arbitrary 
power  law  at  a  specific  wavenumber  is  constant.  (3L  is  found  using  the  following 
relation, 

a(Q)/3LK„-“  -  =  0.  (4.38) 
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Figure  4.3  The  Strehl  ratio  as  a  function  of  a.  PL  is  constant  and  fixed  at  the 
value  specified  in  the  legend.  The  source  aperture  diameter  £)  is  1.0  m. 
The  wavelength  is  1.3  microns. 

Solving  Eq.  (4.38)  for  pL  yields 

PL  =  (4,39) 

cil  ol  j 

This  basis  is  also  appropriate  for  the  ARGUS  data,  where  the  frequency  response  of 
the  anemometer  probes  may  not  corrupt  the  measurement  of  p  in  certain  wavenum¬ 
ber  ranges.  Figure  (4.4)  is  a  plot  of  the  Strehl  ratio  as  a  function  of  the  power  law  a 
for  aim  diameter,  circular  source  aperture.  The  wavenumber  at  which  the  spectra 
are  equal  is  Kq  =  lm“^.  The  wavelength  is  1.3  microns.  One  can  see  from  Figure 
(4.4)  that  the  constraint  forces  PL  to  be  large  when  a  =  3,  thus  po  and  the  Strehl 
ratio  are  small.  As  the  power  law  increases,  the  constraint  forces  PL  to  become 
smaller  until  the  effect  of  pure  tilt  dominates  when  a  4. 
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Figure  4.4  The  Strehl  ratio  as  a  function  of  a.  (31  is  chosen  such  that  Eq.  (4.39) 
is  satisfied  for  the  values  of  C^L  specified  in  the  legend.  The  source 
aperture  diameter  D  is  1.0  m.  The  wavelength  is  1.3  microns.  The 
wavenumber  at  which  the  spectra  are  equal  is  Kg  =  lm“^. 

4.6  Comparing  Strehl  Ratios  in  Arbitrary  Power  Law  Turbulence  Spectra  when  (3L 
is  such  that  the  Power  in  a  Finite  Bandwidth  is  Constant 

In  this  comparison  /?L  is  chosen  such  that  the  power  in  a  finite  bandwidth  is 
constant.  PL  is  found  using  the  following  relation, 

r  d«a(a)^LK-“  -  r  dKa{n/S)ClLK-^^^^  =  0.  (4.40) 

Jki  J 

The  optimal  way  to  accomplish  this  comparison  would  be  to  compute  the  total 
power  in  the  spectrum  instead  of  the  power  in  a  finite  bandwidth.  However,  this  is 
impossible  for  power  law  turbulence  spectra  since  the  total  power  is  infinite.  Figure 
(4.5)  is  a  plot  of  Strehl  ratio  vs.  power  law  for  a  1  m  diameter,  circular  source 
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Figure  4.5  The  Strehl  ratio  as  a  funtion  of  a.  PL  is  chosen  such  that  Eq.  (4.40) 
is  satisfied  for  the  values  of  C^L  specified  in  the  legend.  The  source 
aperture  diameter  D  is  1.0  m.  The  wavelength  is  1.3  microns. 

aperture.  The  bandwidth  in  which  the  spectral  power  is  constrained  to  be  equal  is 
.001  <  K  <  100.  One  can  see  a  similar  effect  on  PL  as  in  the  previous  section.  In 
Figure  (4.5)  the  constraint  forces  pL  to  be  large  when  a  =  3,  thus  po  ^-nd  Strehl  ratio 
are  small.  As  the  power  law  increases,  the  constraint  forces  pL  io  become  smaller 
until  the  effect  of  pure  tilt  again  dominates  when  a  — >  4. 

4.7  Comparing  Strehl  Ratios  in  Arbitrary  Power  Law  Turbulence  Spectra  when  pL 
is  such  that  the  Piston- Removed  Wave  Variance  is  Constant 

The  last  basis  used  to  compare  the  Strehl  ratios  as  a  varies,  is  to  choose  pL 
such  that  the  turbulence  power  spectra  produce  equal  piston-removed  wave  variance. 
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The  wave  variance  is  defined  as 


d+4. 


(4.41) 


where  is  the  log  amplitude  variance  and  (t|  is  the  phase  variance.  Due  to  the 
infinite  power  in  a  general  power  law  spectrum,  the  phase  variance  is  infinite.  Thus, 
Eq.  (4.41)  is  infinite.  However,  if  the  piston  component  of  the  abberations  to 
the  optical  wave  is  excluded,  the  phase  variance  is  finite  and  the  wave  variance  is 
calculable  via  Mellin  transform  techniques.  Thus  we  are  interested  in  the  quantity 


(4.42) 


where  —p  indicates  the  removal  of  the  piston  component. 


4.7.1  Derivation  of  the  piston-removed  wave  variance.  The  spherical  wave 
correlation  functions  for  the  scintillation  and  phase  are  rewritten  from  Eqs.  (2.59 
and  2.60): 


B^{p)  =  47r^A:^  dz  dK«Jo 


k^z[L  —  z) 
2kL 


^^{k,z)  (4.43) 


and 


B*(,)  =  4A^jf"d.fd«y„(^) 


fKpz\  2 

'  ‘  COS 


k^z[L  —  z) 


2kL 


$„(/c,2).  (4.44) 


Since  and  cos2(a;)  +  sin^(a;)  =  1  the  wave  variance  is 


al  =  47r2fc2 


pL  poo 

/  dz  d«:«;$„(K,  z). 
Jo  Jo 


(4.45) 


Equation  (4.45)  is  valid  for  both  plane  and  spherical  waves.  Following  the  approach 
of  Sasiela  [23]  and  Noll  [20],  we  can  remove  the  effects  of  any  Zernike  mode  by 
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introducing  a  filter  function  which  removes  the  effects  of  that  mode.  The  plane  wave 
piston  removed  filter  function  is 


F{k,z)  =  1-4 


Ji(/cT»/2) 


/cD/2  J 


(4.46) 


Multiplying  the  integrand  of  Eq.  (4.46)  by  the  filter  function  gives 


/»L  /»oo  ( 

=  47r^A:^  J  dz  dKK:$„(K:,  2:)  <  1  -  4 


Ji(kD/2) 


kD/2 


(4.47) 


Substituting  the  turbulence  power  spectrum  (Eq.  (2.56))  into  Eq.  (4.47)  yields 


—  ATr‘^k'^a{a)  J  dz(d{z)  J  dKs}  "|l  — 4 


Ji{kDI2) 
kD/2  j 


t2' 


Making  the  change  of  variables 


I  =  —  dx='>''® 

2 

Ox 

K  =  ^  dK  = 


2 

2di 


D 


D  ’ 


(4.48) 


(4.49) 


yields 


£  dzPiz)  £  [1  - 


(4.50) 


If  we  let  s  =  2  —  a,  and  multiply  by  —1,  Eq.  (4.50)  becomes 


= -47r^fc^a(a)2^i9  ''  [  dzP{z)  [  ^  [dx  ^Ji(a;)^  -  ll  .  (4.51) 

JQ  Jo  ^ 


Equation  (4.51)  is  now  in  the  form  of  a  Mellin  transform.  From  the  tables  in  Sasiela 

122|, 


M[x°'h{x)]  H{s  +  a) 
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M[Jl{x)] 


s/2  +  i/  1/2 -s/2 
i/  +  l  —  s/2  1  —  s/2 

-e<^s)<e. 


■2u  <  SR(s)  <  1 


(4.52) 


Applying  the  transforms  in  Eq.  (4.52)  to  Eq.  (4.51)  yields 

al  =  -4Tr^k'^a{a)2W~^-^  [  dzf3{z) 

y/TT  Jo 

X  lime^o  1 


f 

'  +  l  1/2-^ 

1  1  1 

<  r 

9  5-2  1  5-2 

.  2  ^  2  . 

—  — . —  - 

5  +  €  5  — € 

< 

* 

-€  <  3?(s)  <  e.  (4.53) 


In  the  limit  as  e  ^  0  the  poles  at  s  =  0  and  s  —  —e  cancel  and  the  path  of  the 
integration  in  the  complex  plane  moves  to  the  left  one  pole.  Thus,  Eq.  (4.53) 
becomes 


05+1  r)-$ 

alj  =  -47r^fc^a(a) - — E 

w-p  \  J  ^ 


S.  3— s 
2  2 

6—^  4—s 

2  2 


-2  <  ^s)  <  0.  (4.54) 


Since  s  =  2  -  a  the  piston  removed  wave  variance  is 
=  -47r2fc2a(a)  /  P{z)d. 

J  0 


rL  I 

izl  — ^ — r 


-v/tt 


2— g  1— g 
2  2 

4+g  2+g 

2  2 


2  <  a  <  4.  (4.55) 


If  I3{z)  is  constant  over  the  path  L,  the  piston  removed  wave  variance  is 


al_  =  -2^-^TT^I‘^k'^a{a)(3LT 


2— g  1— g 
2  2 

4+a  2+a 

2  2 


2  <  a  <  4.  (4.56) 


4.7-2  Comparing  Strehl  ratios  in  arbitrary  power  law  turbulence  spectra  when 
pL  is  such  that  the  piston  removed  wave  variance  is  constant.  In  this  comparison 
f3L  is  chosen  such  that  the  piston  removed  wave  variance  is  constant.  pL  is  found 
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using  the  following  relation, 


11/3)  =  0.  (4.57) 


Solving  Eq.  (4.57)  for  f3L  yields 


_5  _8  4+a  2+a 

6  6  2  2 

33  16  2:::a 

6  6  2  2  . 

2  <  a  <  4  (4.58) 


Figure  (4.6)  is  a  plot  of  the  Strehl  ratio  as  a  function  of  the  power  law  a  for  a  1  m 
diameter,  circular  source  aperture.  The  wavelength  is  1.3  microns.  One  can  see  in 


pL  =  -2 


-(11/3+Q:) 


q(11/3) 

a(a:) 


(^2_^_pll/3-ar 


a 


— 

P(1  !/3)Z..«i,.4.0xl(r'^ 

— 

KlI/3)Z,.C‘t.2.0xIor‘2 

— 

P(ll/3)t.cji. 4.0x10-'’ 

P(1 1/3)1, .  C'Z,  -  2.0  X  lOr” 

— 

—  • 

P(1 1/3)1.. C'Z,. 4.0x10-“ 

Figure  4.6  The  Strehl  ratio  as  a  function  of  the  power  law.  /3L  is  chosen  such  that 
Eq.  (4.58)  is  satisfied  for  the  values  of  C^L  specified  in  the  legend.  The 
source  aperture  diameter  D  is  1.0  m.  The  wavelength  is  1.3  microns. 
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Fig.  (4.58)  that  the  constraint  on  /3L  now  behaves  differently  depending  the  strength 
of  the  turbulence.  If  the  turbulence  is  weak,  changing  the  power  law  has  little  effect 
under  this  constraint.  If  the  turbulence  is  strong  /3L  can  be  quite  small  at  a  =  3, 
thus  Strehl  ratio  can  approach  1. 


4-19 


V.  Contributions  of  the  Log  Amplitude  and  Phase  Structure 

Functions 

Adaptive  optical  systems  systems  can  only  compensate  for  phase  fluctuations 
in  an  optical  beam.  Therefore,  if  a  distorted  beam  contains  more  phase  fluctuations 
and  less  log  amplitude  fluctuations,  an  adaptive  optic  compensation  system  can  be 
expected  to  perform  a  high  percentage  of  correction. 

In  the  previous  chapters  we  were  concerned  with  the  wave  structure  function 
Dw{p).  It  was  shown  in  Chapter  2  that  the  wave  structure  function  is  actually  the 
sum  of  the  log  amplitude  and  phase  structure  functions.  The  goal  of  this  chapter  is 
to  determine  the  conditions  under  which  either  phase  or  amplitude  perturbations  to 
the  optical  beam  dominate,  as  the  power  law  a  is  varied  from  3  <  a  <  4.  This  will 
be  done  by  plotting  the  relative  contributions  of  the  log  amplitude  D-^[p)  and  phase 
D^{p)  structure  functions  to  the  wave  structure  function  D^[p). 

Rewriting  Eqs.  (2.61)  and  (2.62)  gives 

Dw{p)  =  L>xiP)  +  P<t>ip)  (^-1) 


and 

DM  =  21Bj(0)  -  BM]  (5-2) 


and 

D^{p)  =  2[B^{0)-B^{p)].  (5.3) 

The  log  amplitude  and  phase  correlation  funtions  are  rewritten  from  Eqs.  (2.59  and 
2.60) 


pL  roo 

By^{p)  =  J  dz  J  d«;  kJq 


(Kpz\  .  2 

[-rn 


K^Z 


(L-z) 


2kL 


(5.4) 
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and 


B^{p)  =  47r^A:^  dz  dK  kJq 


cos 


k^z[L  —  z) 


2kL 


^n{i^,z),  (5.5) 


where  Jq  is  a  Bessel  function  of  the  first  kind,  order  0,  k  is  the  optical  wavenumber,  p 
is  the  geometrical  separation  between  points  in  the  source  plane,  and  L  is  the  length 
of  the  optical  path.  By  substituting  Eqs.  (5.4  and  5.5)  into  Eq.  (5.2)  we  get 


pL  poo 

DM  =  dz  dK 


K 


J.(^) 


sm‘ 


'  k^z{L  —  z)' 
2kz 


$„(k,q:,2;)  (5.6) 


and 


pL  poo 

D^{p)  =  S'K^k'^  J  dz  dK 


K 


COS^ 


The  arbitrary  power  law  turbulence  spectrum  from  Eq.  (2.56)  is 


$,,(«,  a,  2)  =  a{a)l3{z)'^  “• 


(5.8) 


Substituting  Eq.  (5.8)  into  Eqs.  (5.6  and  5.7)  yields 


dAp)  = 

87r^A:^a(a)  j  dzf3{z)  J  dK  k^~°‘  1  —  Jo(— ^)  si 


sm 


'  k^z{L  —  2)' 
2kz 


(5.9) 


and 


dAp)  = 


87r^fc^a(a:)  j  dzl3{z)  J  dK  A 


L  L 
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To  derive  some  general  conclusions  about  the  behavior  of  Eqs.  (5.9)  and  (5.10) 
consider  the  change  of  variables 


2:  =  z'L  dz  —  dz'L, 


thus  Eqs.  (5.9)  and  (5.10)  become 


DM  = 

8x^fc^a(tt)  j  dz'f3{z'L)  dK  [1  -  Jo(«p2')] 


and 


D^{p)  =  2  /  / 

87r^k^a{a)  dz'/3{z'L)  dn  /c^"“  [1  -  Joinpz')]  cos^  ^  ^  2fc 


One  can  normalize  the  wavelength  and  path  length  dependencies  by  defining  the 
parameter  Fresnel  ratio  {FR)  which  is  the  ratio  of  the  separation  of  points  in  the 
source  plane  divided  by  the  Fresnel  number  scaled  by  47r 


FR  = 


47rp^  47rp^ 

~\r  ^ 


(5.13) 


Substituting  Eq.  (5.13)  into  Eqs.  (5.11)  and  (5.12)  yields 


D^{p)  =  87r^A:^a(a) 

£  iz'P(z'L)  /;  d.  [1  -  J„(..2')l  (5.14) 

and 
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(5.15) 


D*(p)  =  inVa(a) 

dz'p{z  L)  dK  “  [1  -  JoiKpz  )]  cos  (  - — - j . 

The  problem  is  now  reduced  to  a  double  integration  over  the  turbulence  frequency 
variable  k  and  the  normalized  path  z'.  Eqs.  (5.14)  and  (5.15)  contain  a  product  of 
functions  in  k  which  yield  an  integrand  that  is  highly  oscillatory  in  the  k,  z'  plane 
and  does  not  decay  rapidly.  This  makes  numerical  integration  of  Eqs.  (5.14)  and 
(5.15)  over  both  k  and  z'  quite  dithcult.  Luckily,  both  integrals  can  be  evaluated 
over  K  in  closed  form  using  advanced  Mellin  transform  techniques  according  to  the 
following  procedure: 

1.  Use  the  Mellin  convolution  theorem  to  transform  the  real  integral  into  a  Mellin- 
Barnes  contour  integral  in  complex  plane.  Consider  the  integral 

too 

h{x)  =  /  y~^ho{y)hi{xly)dy.  (5.16) 

JO 

Take  the  Mellin  transform  of  h(x) 

too  too 

M{h{x))  =  /  /  y-'^ho{y)hi{x/y)x^~^dydx.  (5.17) 

V  0  v  0 

By  making  the  change  of  variables  w  =  x/y,  thus  x  =  wy  and  da:  =  ydw  and 
interchanging  the  order  of  integration  Eq.  (5.17)  becomes 

poo  too 

M{h{x))  =  /  y  ho{y)dy  /  {wyy~^hi{w)ydw.  (5.18) 

1/  0  v  0 

Simplifying  Eq.  (5.18)  yields  an  integral  of  the  Mellin  transforms  of  hg  and  hi 

poo  poo 

M(h(x))  =  {wY-%iw)dw  =  (5.19) 
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The  inverse  Mellin  transform  is  now  the  contour  integral 


h{x)  Ho{s)H^{s)x-^ds,  (5.20) 

Ji'KJ  Jc 

where  Ho{s)  and  i?i(s)  are  the  individual  Mellin  transforms  of  ho{y)  and 
hi{x/y).  Since  the  Mellin  transforms  of  ho{y)  and  hi{xfy)  in  Eq.  (5.20)  can  be 
expressed  as  ratios  of  gamma  functions,  the  original  integral  h{x)  (Eq.  5.16) 
can  be  expressed  as  a  Mellin-Barnes  integral  of  the  following  form: 

,  /  X  _  1  f,  -s  n^iElai  +  ais]U.f^iT[bj  —  ^js]  Cpioil 

“  2^j  I  nt,ricj  +  7i«]ns=,r|d„  +  «„*]'  '  ’  ' 

2.  Close  the  path  of  the  contour  integral  at  infinity,  subject  to  the  conditions 
specified  in  [22]  and  perform  the  contour  integration  using  Cauchy’s  residue 
theorem.  Application  of  the  residue  theorem  yields  an  infinite  sum  of  the 
residues  at  the  poles  of  the  gamma  functions  in  Eq.  (5.21) 


h{x)  =  '^x  ^”i?es(s„).  (5.22) 

n 

The  poles  occur  only  when  s  —  — (n  +  ai)fai  or  s  =  — (n  +  bj)lf3j. 

3.  Use  the  definition  of  a  generalized  hypergeometric  function  to  convert  the  in¬ 
finite  sum  found  in  Step  2  into  a  sum  of  generalized  hypergeometric  functions. 
The  definition  of  a  generalized  hypergeometric  function  is  given  by  [22] 


p 


OO 


Fg{a,b,z)  = 

k=0 


^  (®l)fc(®2)fc  •  •  •  (®jj)fc 

Id  {bi)k{b2)k...{bi)k' 


(5.23) 


where  (a)  =  ai,  02, . . . ,  Op,  (5)  =  5i,52,-.-,5g  and  {a)k  is  the  Pochhammer 
symbol  defined  as 


(a)fe  = 


r[a  +  k] 
r[a] 


(5.24) 
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Although  the  procedure  outlined  above  seems  quite  difficult  to  the  uninitiated, 
the  method  reduces  to  a  series  of  straightforward  and  logical,  but  somewhat  tedious 
steps.  The  details  are  shown  in  Sasiela  [22]. 

Fortunately  in  this  case  there  is  a  short  cut,  the  k  integration  in  Eqs.  (5.14)  and 
(5.15)  can  be  accomplished  using  Mathematica’s  [29]  symbolic  integration  routines. 
Unfortunately  Mathematica  does  not  always  produce  symbolic  output  in  the  most 
convienient  form,  either  for  physical  interpretation  or  numerical  condition.  In  the 
case  of  Eqs.  (5.14  and  5.15)  the  Mathematica  solution  is  sufficiently  well  behaved  to 
allow  numerical  integration  over  z'  for  the  ranges  of  the  parameters  of  interest  {FR 
and  O').  The  results  of  the  k  space  integrations  performed  using  Mathematica  are 
contained  in  unreduced  form  in  Appendix  A. 

Once  the  integration  over  k  is  performed,  the  sum  of  generalized  hypergeomet¬ 
ric  functions  created  in  Step  3  can  be  numerically  integrated  over  the  normalized 
path  to  get  the  log  amplitude  and  phase  structure  functions. 

Since  we  are  concerned  with  the  relative  contributions  of  the  log  amplitude 
and  phase  structure  functions  to  the  wave  structure  function,  we  can  examine  the 
ratios  D^{p)lD^{p)  and  D^{p)lD^{p)  which  are  given  by 

dm  _  /o^  dz'Pjz'L)  /o°°  dK  k}--^  [1  -  JojKpz')]  sin^ 

DM  ~  So  dz'/^(z'L}  /o°°  dK  [1  -  Mkpz')J 

and 

DM  _  So  dz'/3(z'L)  Jo°°  dK  [1  -  Jo(Kpz')]  cos^ 

DM  ~  So  dz'/5(z'L)  /o°°  dK  [1  -  M/^pz')] 

Close  examination  of  the  equations  for  D^{p)  and  D^{p)  in  Appendix  A  shows  both 
contain  a  dependence  on  The  wave  structure  function,  Eq.  (2.70)  derived 

in  Chapter  2,  also  contains  a  p°‘~^  dependence.  Therefore  ratios  D^{p) / D^{p)  and 
D^{p)lDw{p)  do  not  depend  on  p  directly.  However,  the  ratios  D^{p)fD^{p)  and 


(5.26) 


(5.25) 
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D^{p)l Dw{p)  do  depend  on  p  indirectly  through  the  definition  of  Fresnel  ratio  in  Eq. 
(5.13). 

Figure  (5.1)  contains  a  plot  of  the  ratios  D^{p)  and  D^{p)  as  a  function  of 
the  power  law  a  for  various  Fresnel  ratio’s  FR.  The  upper  family  of  curves  are 
the  phase  contributions  D^{p)  to  the  wave  structure  function  Dw{p)  and  the  lower 
family  of  curves  are  the  log  amplitude  contributions  to  the  wave  structure  function 
D-u,{p).  The  sum  of  the  log  amplitude  and  phase  contributions  equals  1.  Two  key 


Figure  5.1  The  relative  contributions  of  the  log  amplitude  and  phase  structure 
functions  to  the  wave  structure  funtion  as  a  function  of  the  power  law 
a.  The  Fresnel  ratio  FR  is  fixed  at  the  value  specified  in  the  legend. 

observations  can  be  made  from  this  plot.  Phase  perturbations  increase  as  Fresnel 
ratio  increases.  This  is  consistent  with  the  definition  of  Fresnel  ratio.  For  a  fixed 
separation  p  in  the  source  plane,  the  phase  effects  get  large  as  the  Fresnel  number 
decreases.  In  other  words  phase  effects  increase  when  the  optical  path  is  short.  The 
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second  and  most  important  observation  is  that  phase  effects  dominate  as  the  power 
law  increases  to  o;  =  4.  This  is  consistent  with  the  work  of  Heidbreder  [12]  who 
showed  that  a  wave  structure  function  which  had  a  dependence  corresponded  to 
a  pure  random  tilt  in  the  aperture  of  an  incoherent  imaging  system. 


VI.  Results,  Conclusions,  and  Recommendations 

This  chapter  summarizes  the  major  results  of  this  research.  Following  the 
summary,  the  implications  for  adaptive  optical  beamforming  systems  are  discussed 
and  recommendations  for  future  research  are  made. 

6. 1  Summmary  of  Results 

In  Chapter  2,  Section  2.2,  the  fundamental  relationship  between  the  index 
structure  function  and  the  index  power  spectrum  were  examined.  Assuming  the 
turbulence  is  isotropic  and  the  index  structure  function  obeys  an  arbitrary  power 
law 

D.(p)  =  I3(z)p"~^  (6-1) 

the  corresponding  index  power  spectrum  is 

#„(«;)  =  o(a!)/?(2)K~“,  (6.2) 

where  the  function  a{a)  maintains  consistency  between  the  index  structure  function 
and  the  index  power  spectrum  and  is  given  by 

a{a)  =  2““®(a^  —  5o!  +  3  <  a  <  5.  (6.3) 

^  L~. 

Inner  and  outer  scale  effects  were  neglected. 

In  Sections  2.3  and  2.4  the  wave  structure  function  was  derived  using  the  log 
amplitude  and  phase  correlation  functions  which  can  be  derived  from  the  method  of 
small  perturbations.  The  wave  structure  function  is  given  by 

DM  =  Cl  ,  (6.4) 
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where  Ci  is  given  by 


Cl  =  2 


a-2 


L«-2J 


a-2 

2 


(6.5) 


and  a  new  generalized  definition  for  the  coherence  diameter  is 


Po  = 


c.(a-l)r  [f] 

-(2)4-“7r2A:2a(a)L2-“r 

3  <  a  <  4.  (6.6) 


The  important  conclusions  to  draw  from  Sections  2.3  and  2.4  are,  if  inner  and 
outer  scale  effects  are  neglected,  only  power  laws  from  3  <  a  <  4  are  allowed  to 
exist.  The  physical  interpretations  for  the  mathematical  results  are: 


1.  When  the  power  law  o:  — ^  3,  the  index  structure  function,  Eq.  (6.1),  has  no 
dependence  on  p,  therefore  no  turbulence  exists. 

2.  When  the  power  law  a  — >  4,  the  wave  structure  function,  Eq.  (6.4),  has  a 
dependence  and  the  resulting  abberation  to  the  optical  wave  is  a  random  tilt 
in  the  source  aperture  which  causes  the  beam  to  wander,  thereby  dramatically 
increasing  the  long  exposure  spot  size. 

In  Section  2.5  the  variance  of  the  log  amplitude  fluctuations  was  found  to  be 


a 


2 

X 


4ci7r*2 


a 


(6.7) 


where  FN  is  the  Fresnel  number  defined  as 


FN  =  a/AL.  (6.8) 

The  log  amplitude  variance  is  a  maximum  as  a  — >  3.  The  value  of  the  maximum  log 
amplitude  variance  is  determined  by  the  ratio  of  the  Fresnel  number  to  the  coherence 
diameter  [FNj Po).  The  log  amplitude  variance  decreases  greatly  as  a  — >  4  and  phase 
effects  dominate. 
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In  Chapter  3,  the  mutual  coherence  function  of  an  arbitrary  wave  propagating 
through  turbulence  which  obeys  an  arbitrary  power  spectrum  was  derived.  The 
result  is  a  4-dimensional  integral  of  the  form 


MCF{xi,yi,X2,y2)  = 

H{xuyu  6,  ni,  ^2,  y2, 6,  z) 

X  exp  +  rji)  -  (^1  +  'nl)  +  22:26  -  2x16  +  22/2^2  -  2yir7i) 

xC^a(6,  ^2)d6dmd6d^72, 


(6.9) 


where  H(xx,yi^  6?  'Qu  2^2, 1/2, 6)  ^2,  z)  is  the  second  moment  of  the  Greens  function  for 
propagation  of  a  spherical  wave  in  aturbulent  medium.  H(xi,yi,^i,r]i,X2, 2/2, 6)  V2,  z) 
is  shown  in  [17]  to  be 


H{xi,yi,^ur)i,X2,y2,^2,V2,z)  =  exp 


(6.10) 


From  the  general  form  of  Eq.  (6.9)  many  parameters  of  interest  can  be  calculated, 
including  intensity  profiles  and  Strehl  ratios  for  arbitrary  apertures. 

The  Strehl  ratio  of  a  constant  amplitude  beam  is  derived  in  Chapter  4.  The 
MCF  derived  in  Chapter  3  was  reduced  to  a  single  numerical  integral  of  the  form 


STZ  =  —  [  -(cos  ^(u)  —  uy/\  —  V?)  exp 
TT  20  2 


1  J 

(D 

— u 

2  ' 

\Po  ) 

udu. 


(6.11) 


Several  comparisons  of  the  Strehl  ratio  as  a  function  of  the  power  law  a  were  made 
under  various  constraints.  It  was  necessary  to  use  a  variety  of  constraints  to  compare 
the  Strehl  ratio  at  power  law  ai  to  Strehl  ratio  at  another  power  law  0:2  because 
changing  a  changes  the  total  power  in  the  turbulence  spectrum,  therefore  biasing 
the  results  one  way  or  another.  In  general,  the  resulting  trends  of  the  Strehl  ratio 
depend  on  the  constraint  placed  on  the  integrated  turbulence  strength  /3L  as  the 
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power  law  a  was  varied.  The  consistent  trend  noticed  in  the  behavior  of  the  Strehl 
ratio  was  that  the  Strehl  ratio  approached  zero  as  o:  — >  4,  this  behavior  is  consistent 
with  the  physical  model  of  a  pure  tilt  aberration  in  the  source  aperture. 

In  Section  4.2  the  asymptotic  behavior  of  Fried’s  UjUmax  resolution  metric 
was  examined  for  a  general  power  law.  The  design  of  this  metric  is  such  that  when 
the  power  law  ot  is  varied  from  3  <  ct  <  4,  it  is  possible  to  define  the  generalized 
coherence  diameter  po  and  the  constant  Ci,  such  that  the  interpretation  of  po  is  the 
same  as  Fried’s  coherence  diameter  for  Kolmogorov  turbulence,  to- 

Finally,  in  Chapter  5  the  relative  contributions  of  the  log  amplitude  and  phase 
structure  functions  to  the  wave  structure  functions  were  calculated.  This  analysis 
showed  conclusively  that 

1.  As  the  power  law  o:  ^  3  the  relative  contributions  of  the  log  amplitude  and 
phase  fluctuations  are  soley  controlled  by  the  Fresnel  ratio  of  the  system. 

2.  As  the  power  law  a  A  phase  effects  dominate  the  wave  structure  function 
regardless  of  Fresnel  ratio. 

It  is  postulated  from  the  results  of  Chapters  4  and  5  that  if  the  perturbations  of  the 
wavefront  in  the  source  plane  is  decomposed  into  a  series  of  Zernike  polynomials, 
the  energy  in  the  series  shifts  toward  the  low  order  modes  when  the  power  law  a 
increases  towards  4. 

6.2  Implications  for  Adaptive  Optical  Beamforming  Systems 

An  adaptive  optical  beamforming  system  corrects  for  atmospheric  effects  by 
predistorting  the  phase  of  the  wavefront  to  be  transmitted  with  the  conjugate  of 
the  phase  perturbations  induced  in  the  atmosphere.  Thus,  atmospheric  conditions 
which  yield  only  phase  perturbations  are  more  conducive  to  the  operation  of  adap¬ 
tive  optical  systems.  This  means  turbulence  power  spectra  which  have  power  laws 
approaching  4  are  more  conducive  to  the  operation  of  adaptive  optical  systems.  For 
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equal  turbulence  strength,  the  ABL  system  can  expect  better  performance  in  turbu¬ 
lence  whose  power  law  tends  toward  a  =  4.  The  converse  is  also  true,  regardless  of 
the  Fresnel  ratio  of  the  system,  the  ABL  system  can  expect  worse  performance  as 
O!  — >  3. 

Additionally,  if  the  limit  on  the  method  of  small  perturbations  (a^  <  0.3  -  0.5) 
is  assumed  to  hold  regardless  of  the  power  law,  the  method  of  small  perturbations 
can  be  used  at  longer  path  lengths  and  greater  turbulence  strengths  than  previously 
determined  if  o:  — >  4. 


6.3  Recommendations  for  Future  Research 

In  this  section  several  extensions  of  this  research  are  proposed. 

1.  Determine  the  frequency  response  of  the  ARGUS  anemometer  probes.  This  is 
the  highest  priority,  because  without  an  accurate  calibration  of  the  frequency 
response  the  ARGUS  anemometer  data  is  inconclusive.  With  the  current  data 
it  is  impossible  to  separate  the  effects  of  inner  and  outer  scales  of  turbulence 
from  the  power  law  of  the  inertial  subrange  and  from  the  probe  response.  Once 
the  calibration  has  been  accomplished,  these  effects  are  separable  and  one 
could  determine  the  directions  of  further  investigations  into  non-Kolmogorov 
propagation  theory. 

2.  Perform  a  modal  decomposition  of  the  propagating  wavefront  using  Zernike 
polynomials.  By  estimating  the  average  energy  associated  with  each  polyno¬ 
mial  in  the  series,  the  postulate  in  Section  6.2  can  be  proved  or  disproved. 

3.  The  analysis  presented  in  this  thesis  neglected  the  effects  of  an  inner  and  outer 
scale.  A  generalized  power  spectrum  incorporating  the  effects  of  an  inner  and 
outer  scale  can  be  defined  as 


a,  z)  =  a'{a)/3(z){K^  -1-  kI)  exp 


(6.12) 
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where  and  Ko  are  the  inner  and  outer  scales,  respectively.  Note:  a'{a)  in 
Eq.  (6.12)  is  not  the  a{a)  defined  in  this  thesis. 
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Appendix  A.  Integration  of  the  Log  Amplitude  and  Phase  Structure 

Functions 

In  Chapter  5,  Mathematica  was  used  to  evaluate  the  following  integrals.  The 
integrals  over  the  turbulence  wavenumber  k  were  evaluated  in  closed  form.  The 
integrals  over  the  path  z'  were  evaluated  numerically.  Specifically,  from  Eqs.  (5.14 
and  5.15)  we  have 

Pxip)  ~  8Tr'^k'^a(a) 

X  ^  dz'f3{z'L)  d/c  [1  -  Jq{kpz')]  sin^ 

and 


D^{p)  =  87r^A:^a(a) 


X  J  dz'l3{z'L)  j  d/c  “  [1  -  Jq{kpz')] 


cos 


[Kpfzfl  —  z') 

Jr 


(A.2) 


Integrating  Eq.  (A.l)  symbolically  over  k  gives 

Dy^{p)  —  SJk^a{a)  J  dz'P{z'L) 

^  _ ^FR,^ 

if  csc(^)pF,({l-f  ,|-f 


2  ) 


128^f  2f  (-l+f)p2(l-^)"f^*  r(-l+f) 

^  a  {p^  ^  csc(^)pF^({l-f  ,|-f  ) 


,2  FR^  4  p/  1  I 
3  a 


5125!^^  2f  (-H-f)p2(l-A 

3if  7r(p^z^)f  secj^) 

128  2f  (i-f  )  (l-f  )  (-l  +  f  )  FR^  p^  z^  r(-2+f ) 

^  lfa.(p2,2)f^|-f  3ec(^) 

256  ^  2f  (i-f)  (l-f)  (-l+f)  FR?  r(-2+f ) 

64  2f  (i-f)  (l-f)  (-l+f)  FR^  p2  z3  r(-2+f 


+ 


+ 


(A.3) 
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128  2f  (i-f  )  (l-f  )  (-l+f  )  FR^  r(-2+f ) 


— I  KR^ 

o  1  ^  /  .2  .2  \  2  FR^  4  m  /(I  ori  aifllit  \ 

3z^Tr{pz)  vFM2  4’^  4  J>t2’2’^>’  256  (1-z)2  ^ 

128  5ij»  2f  (l-f  )  (l-f)  (-l+f)  FR^  P^  r(-2+f ) 


)  sec(^) 


256^^  2f  (l-f)  (l-f)  (-1 +f  )  FR^  p^  z^  r(-2+f ) 

where  pFq{)  is  the  generalized  hypergeometric  function.  Integrating  Eq.  (A. 2) 
symbolically  over  k  gives 


D^{p)  =  SF^k‘^a{a)  J  dz'l3{z'L) 


«  «  1  _L  Ck: 

2  ,2  -7+¥  , 


.  (1-^)2  ^2-?+f  . 


2(2-a)(l-f)(l-2)^  r(-l+f)  (2-a)  (l-f)  (1-^)2  r(-l+f) 


2  .2  -Hf 


2(2-«)  l-f  (i-^)^r(-i+f) 


r(|-f)  I 

/  ■•  .  rv  \  I  / 


-f)7r(p=^z^)^  csc(^) 


(-l+f)p2^2r(a)2 


(A.4) 
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